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ABSTRACT 

An  experimental  investigation  was  conducted  to  examine 
fuel  regression  rate  control  methods  other  than  variable 
bypass  air  flow  rates  in  the  solid  fuel  ramjet.  Air  and 
oxygen  injection  at  various  axial  locations  within  the  fuel 
grain  were  examined  as  well  as  air,  oxygen  and  ethylene 
injection  through  the  step  face.  One  inlet  swirl  design  was 
also  tested.  Secondary  gas  injection  was  found  to  be 
inadequate  for  regression  rate  control.  A  small  amount  of 
inlet  swirl  resulted  in  a  significant  increase  in  fuel 
regression  rate,  indicating  that  variable  inlet  swirl  may  be 
a  viable  technique  for  providing  in-flight  fuel  flow  rate 
modulation  in  the  solid  fuel  ramjet. 
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I.  INTRODUCTION 

Ramjets  operate  with  much  higher  specific  impulse  than 
rockets  since  ramjets  use  inlet  air  as  a  source  of  oxygen. 
Self-sufficient  rockets,  on  the  other  hand,  must  carry  their 
own  oxidizer  and  bear  the  consequent  penalty.  Accordingly, 
although  rockets  must  be  chosen  for  propulsion  outside  of 
the  atmosphere,  and  solid  propellent  rockets  are 
unchallenged  for  short-  range  tactical  applications,  ramjets 
can  generally  outperform  rockets  in  the  medium  and  long 
range  tactical  environments. 

Because  the  ramjet  depends  only  on  its  forward  motion  at 
supersonic  speeds  to  effectively  compress  intake  air,  the 
engine,  in  principle,  can  employ  very  few,  if  any,  moving 
parts.  It  is  therefore  capable  of  simplicity,  lightness  of 
construction,  and  high  flight  speed  not  possible  in  other 
air-breathing  engines.  These  features,  plus  the  high  thermal 
efficiency  it  can  achieve,  make  the  ramjet  a  particulary 
attractive  choice  for  propelling  vehicles  at  supersonic 
speeds . 

One  of  the  significant  difference  between  rockets  and 
ramjets  is  thrust  at  zero  speed.  Rockets  can  deliver  thrust 
at  any  speed,  whereas  a  ramjet  requires  an  auxilary  boost 
system  to  accelerate  it  to  its  supersonic  operating  regime 
so  that  its  forward  motion  can  compress  the  inlet  air.  To 
operate  at  practical  efficiency  a  ramjet  must  be  moving  at 
about  a  Mach  number  of  1.5  or  greater  so  that  the  margin  of 
thrust  over  drag  will  be  satisfactory. 

The  solid  fuel  ramjet (SFRJ)  employs  solid  fuel  for  the 
combustor  walls.  Its  distinguishing  feature  is  the  absence 
of  fuel  tankage,  fuel  delivery,  and  fuel  control  systems. 
The  solid  fuel  ramjet  offers  this  simplicity  while  providing 
excellent  density  impulse  and  combustion  efficiency. 
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Fig  1.1   shows  a   schematic  of   a  SFRJ   combustor.  One 

inherent  problem   with  the   SFRJ  is   the  dependence   of  fuel 

regression  on  the  air  mass  flux  through  the  fuel  grain.  The 
fuel  regression  rate  generally  behaves  according  to 


,PJ  (1.1) 


where  f     =  fuel  regression  rate 
Mair  =  air  flow  rate 
Ap    =  port  area 

f(Tair,P)=   a    weaker   function    of   inlet    air 
temperature  and  pressure 

x   is  between  . 3  and  . 6 

If  the  air  flow  rate  increases  or  decreases,  then  the 
fuel  flow  rate  varies  in  the  correct  direction,  but  not 
enough  to  maintain  design  fuel-air  ratio.  This  can 
significantly  affect  propulsive  thrust. 

The  range  of  application  and  the  performance  of  the  SFRJ 
could  be  significantly  improved  if  fuel  regression  rate 
could  be  controlled  in  some  manner.  Perhaps  the  most 
obvious  method  for  fuel  regression  rate  control  is  to  use 
variable  bypass  air  flow. 

Bypass  designs  are  often  desired  in  order  to  increase 
fuel  loading,  and  can  also  improve  combustion  efficiency.  A 
valve  in  the  bypass  line  could  be  used  to  vary  the  bypass 
ratio,  thus  producing  changes  in  the  fuel  flow  rate 
according  to  equation  1.1.  (since  only  the  fuel  port  air 
mass  flux  affects  the  regression  rate).  However,  combustion 
efficiency  may  vary  significantly  with  the  amount  of  bypass 
air  and  a  control  valve  must  be  used  in  the  inlet  air 
ducting.   In  addition,   flow  coupling   may  lead  to  undesired 
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oscillations.  Therefore,  other  fuel  regression  rate  control 
techniques  should  be  persued. 

One  possible  alternative  technique  is  the  use  of 
variable  swirl  at  the  air  inlet.  Increased  swirl  may 
increase  the  regression  rate.  This  technique  would  require  a 
vane  control  device.  Another  possibility  is  the  use  of 
secondary  injection  of  air,  oxygen  or  gaseous  fuel  into  one 
or  more  locations  within  the  combustor.  This  may  lead  to 
both  increased  fuel  regression  rate  and  combustion 
efficiency.  If  oxygen  were  used  it  would  require  an 
auxiliary  supply  system.  Therefore,  it  would  be  a  practical 
possibility  (to  maintain  simplicity)  for  a  one-time 
augmentation  in  thrust(at  take-over  from  booster,  at 
terminal  maneuver,  etc). 

In  this  investigation  tests  were  conducted  to  examine 
the  effects  of  secondary  gaseous  injection  and  inlet  air 
swirl  on  regression  rate  and  combustion  efficiency. 

A  new  ramjet  thrust  stand  was  installed  in  the 
combustion  laboratory  and  required  calibration  and 
certification.  A  second  part  of  this  investigation  was  to 
conduct  tests  using  high  temperature  inlet  air  and  then  to 
ensure  that  combustion  efficiencies  based  on  thrust  and 
based  on  measured  chamber  pressure  were  in  agreement. 
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II.  METHOD  OF  INVESTIGATION 

A.   TEST  CONDITITIONS  AND  METHODS 

The  experimental  investigation  was  begun  by  selecting 
several  configurations  for  regression  rate  control  through 
the  utilization  of  secondary  injection.  It  was  decided  to 
examine  secondary  injection  at  three  locations  through  the 
fuel  grain  wall:  in  the  recirculation  zone(flame  holder 
area),  just  downstream  of  flow  reattachment  and  further 
downstream  within  the  region  of  the  developing  boundary 
layer  (in  which  is  located  a  turbulent  diffusion  flame).  In 
addition,  gaseous  injection  on  the  inlet  step  face  and  inlet 
air  swirl  were  examined. 

Polymethy lmethacry late (PMM)  was  used  as  the  fuel  and  had 
a  nominal  length  of  12  inches  and  an  internal  diameter  of 
1.5  inches.  A  sudden  expansion  inlet  was  used  which  provided 
an  inlet  step  height  of  0.49  inches. 

Flow  reattachment  normally  occurs  between  7  and  8  step 
heights (i.e,  3.5  -  4  inches).  Thus  the  three  injection 
location  chosen  where  2,  5,  and  8  inches  from  the  head 
end(See  Figure  2.1). 

The  injection  ring  was  designed  to  provide  radial 
injection  through  eight  equally  spaced,  0.0625  inches 
diameter  holes (See  Figure  2.2  and  Figure  2.3). 

Gaseous  injection  was  also  used  through  the  inlet  step 
face  as  shown  in  Figure  2.4  and  Figure  2.5. 

Swirl  was  provided  to  the  inlet  air  flow  through  the  use 
of  a  swirl  element  inlet  inj ector (Figure  2.6). 

Initial  screening  tests  of  short  duration  were  made  in 
which  no  injection  or  swirl  was  used,  followed  by  tests  with 
various   amounts   of   gaseous   injection   or   swirl.    Those 
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conditions  found  to  significantly  affect  chamber  pressure(or 
thrust)  were  further  evaluated  using  full-length  burn  times 
of  approximately  30  seconds  (in  order  to  accurately 
determine  regression  rate  and  combustion  efficiency). 
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DATA  COLLLECTION  METHOD 


The  data  collected  during  the  regression  rate  control 
tests  consisted  of  air  flow  rates,  oxygen,  air  or  fuel  flow 
rates  for  the  injection,  motor  head-end  and  chamber 
pressures,  weight  changes  of  the  fuel  grain,  pre-ignition 
air  flow  time,  ignition  time,  burning  time,  purge  time, 
thrust  and  ignition  gas  flow  rate.  All  pressures  and  the 
thrust  were  recorded  both  digitally  and  analog. 


Figure  2.1    Fuel  Configurations  for  Side  Wall  Injection 
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Figure  2.3    Injection  Ring  for  Side  Wall  Injection 


Figure  2.4    Step  Insert  For  Gaseous  Face  Injection 
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Figure  2.6    Swirl  Element  Inlet  Injector. 

Temperatures  were  only  recorded  digitally.  Average 
values  of  the  thrust  and  chamber  pressure  were  determined 
using  a  compensating  polar  planimeter  on  the  analog  record. 

The  mass  flow  rate  of  air  was  obtained  by  using 
sonically  choked  nozzles.  Average  fuel  mass  flow  rate  was 
calculated  by  using  the  weight  loss  during  the  burn  time 


Mf  = 


a.  W 


tb 


(2.1) 


where  ^w  •=  weight  change  during  the  run 
tb=  burn  time 
The   average  internal   diameter  of   the   fuel  grain   was 
measured  prior  to  the  run.     The  final  average  diameter  was 
determined  based  on  weight  loss  and  fuel  length  by  using 
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H^fr-  +  51*  (2-2) 

where  AW  =  weight  change 

L   =  length  of  the  fuel  grain 
$      =    density  of  the  fuel  grain 
D]J  =   average  intial   port  diameter   of  the 
PMM  grain 


The  average  fuel  regression  rate  was  then  computed  using 


Pf  -  Pi 

"i  z  tb 


rcwo  "   — (2-3) 


The  mixture  ratio,  chamber  pressure  and  motor  air  inlet 
temperature  were  used  as  input  into  the  Naval  Weapons 
Center(NWC)  China  Lake,  Ca.,  Propellant  Evaluation 
Program(PEPCODE )  computer  program  to  obtain  the  theoretical 
adiabatic  combustion  temperature  and  the  combustion  gas 
properties ( T and  R) .  This  temperature  was  used  to  calculate 
temperature  rise  combustion  efficiencies  based  on  thrust  and 
based  on  nozzle  stagnation  pressure,  where 


% 


T 


Tte*P I* (2.4) 


:ure 


and    T^.^   =  combustor  stagnation  temperati 

T-^j,    =  theoretical  combustor  temperature 
Ta     =  inlet  air  temperature 


Experimental  values  of  combustor  stagnation  temperatures 
were  calculated  in  two  ways ,   one   based  on  pressure  and  one 
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based  on  thrust.  The  one  dimensional  continuity  equation 
(m=  9  A  V)  can  be  expressed  in  terms  of  the  chamber 
stagnation  properties  for  a  choked  converging  nozzle 


•   =  _^^r    g.Y_J.  (2.5) 


where  mi  =  total  flow  rate 

Pt   =   Pc,    chamber   pressure(for   low   Mach 
number) 

Athpfr^  effective  nozzle  throat  area 
R   =  gas  constant 


The  characteristic  exhaust  velocity,  c    is  defined  as 


f=     Pc   Ath  %    __      F  3c 

TT\t  Cp    lYIf 


(2.6) 


where  Cp  =  thrust  coefficient 
For  a  sonically  choked  exhaust  nozzle, 


C*  =  /-|c  (-3—-)1  '   R  Tt  (2.7) 


Using  pre-test  air  flow  through  the  motor,  equations 
(2.6)and  (2.7)  can  be  used  to  determine  the  effective 
exhaust  nozzle  throat  diameter (Dtheff) 
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Dt 


Nf 


C  air   THgi  r 


(2.8) 


Then   temperature- rise   comustion   efficiency   based   on 
nozzle  stagnation  pressure  can  be  calculated  using  c* 


% 


^  exp 2:  a- 

r*  *         r*z 


(2.9) 


where  Ca  =  characterit ic   exhaust  velocity  for  air 


flow  before  ignition 
Cexp   : 


experimental   characteristic   exhaust 


velocity  based  on  Pc,  the  average  combustion  pressure 

C  +.l  is  obtained  from  equat  ion(  2  .  7  )  using  ~Y  , 
R  and  Tt  from  PEPCODE  with  inputs  of  Pc,  Ta,  ma,r,  mi,  m0iand 
mCxH4  •  ^e  latter  two  flow  rates  are  used  in  the  vitiated 
air  heater. 


The  thrust  equation  is 


F    =      ilit    Ue    +     (Pe    -    Po)    Ae 


(2.10) 


*  • 


where      m*  =n^  +    m^+    m0  +    ™UW<± 


A  thrust  coefficient ,  Cp,  can  be  defined  such  that 


F  =  Cp  Pt  Ath 


(2.11) 


^i-l 


Yi 


ex 


2^ 


)(  )     h-(-4J-) 
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J 


Pc   ** 


(2.12) 


where  Pe  =  exit  pressure  of  the  nozzle 

P©  =   ambient  pressure.  Then  for   a  choked 
converging  nozzle,  this  equation  can  be  simplified  to 


where  R  and  j  are  determined  using  PEPCODE 

Then  efficiency  based  on  thrust  is  calculated  from 
*  z    r  *2 

Ce*f> W_  (2    14) 


7. 


«  v.  th      -     Co. 


ignition) 
record 


where    Cexp         =    F    gc  /  (Wj    Cpexp) 


,*- 


Ca  =  F  /  (Mpt*  Cp)  (determined  before 
F  =  average  thrust  from  analog 
Mp*     =  pre-ignition  flow  rate 
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III.  DESCRIPTION  OF  APPARATUS 

A.   RAMJET  MOTOR 

The  ramjet  motor  assembly  used  in  this  experiments  was 
that  used  previously  at  NPS  [Ref.  1,  2]. 

Figure  3.1  shows  a  schematic  drawing  of  the  ramjet, 
illustrating  the  main  sections.  These  are,  the  head-end 
assembly,  step  insert  section,  fuel  grain,  aft  mixing 
chamber  and  exhaust  nozzle.  For  the  tests  with  secondary 
injection,  the  injection  ring  was  used. 

The  head-end  assemly  contains  a  central  opening  for  the 
introduction  of  the  primary  air-flow,  and  ports  for 
introduction  of  the  ignition  fuel  and  the  igniter  torch. 

The  fuel  wall  injection  is  shown  in  Figure  2.1,  Figure 
2.2  and  Figure  2.3.  Injection  velocity  varied  between 
approximately  30  and  200  ft/sec  as  the  injection  mass  flow 
rate  was  increased  from  1%  to  5  %  of  the  inlet  air  flow.  A 
nominal  one  dimensional  flow  port  velocity  was  approximately 
300  ft/sec. 

For  the  face  injection  tests,  air,  oxygen  or  gaseous 
fuel  was  introduced  using  a  differently  designed  step 
insert (Figure  2.4  and  Figure  2.5).  The  injection  was 
provided  through  eight  equally  spaced,  0.047  diameter  holes. 
For  an  injection  gas  mass  flow  rate  equal  to  1  %  of  the 
inlet  air  flow  rate  the  injection  velocity  was  approximately 
65  ft/sec. 

To  provide  a  swirl  at  the  air  inlet,  a  specially 
designed  insert  was  used(Figure  2.6).  This  device  was  a 
tube-in-hole  type  injector.  Approximately  43  %  of  the  air 
flow  passed  through  the  outer  annulus  to  maintain  the 
recirculation  region/flame  holder.  The  other  57  %  of  the  air 
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flow  was  passed  through  a  central  swirl  element  which 
extended  0.5  inches  into  the  combustor.  The  swirl  was 
induced  by  machining  six  flutes  with  a  twist  of  one  turn  in 
five  inches .  The  latter  resulted  in  an  inlet  swirl  angle  of 
approximately  5  degrees  to  the  motor  centerline. 

During  the  normal  ignition  sequence,  the  air  flow  rate 
was  set  at  the  desired  value  before  the  ignition  switch  was 
activated.  When  the  ignition  circuit  was  activated,  ethylene 
gas  was  introduced  upstream  of  the  fuel  grain  along  with 
primary  air  flow.  This  mixture  was  ignited  with  an 
oxygen- ethylene  torch  which  issued  from  the  face  of  the  step 
inlet.  Normally,  two  to  three  seconds  of  ignition  time  was 
required  for  PMM  combustion  to  sustain. 

The  step  insert  section  was  provisioned  for  variations 
in  inlet  diameters  to  make  the  sudden  expansion.  The  fuel 
grain  itself,  when  mounted  in  the  motor,  became  the  mid-body 
of  the  ramjet.  The  fuel  grain  (and  injection  ring,  if  used) 
mounted  between  the  head-end  assembly  and  the  aft  mixing 
chamber.  The  aft-mixing  chamber  had  a  length-to-diameter 
ratio  of  2.9.  The  entire  motor  was  held  together  by  four 
threaded  rods  and  nuts.  The  ramjet  motor  was  then  mounted 
on  the  thrust  stand. 

B.   AIR  SUPPLY  SYSTEM 

A  schematic  of  the  ramjet  air  supply  system  is  shown  in 
figure  3.2.  From  the  air  tanks,  the  air  flows  through  a 
pressure  regulator,  a  sonic  choke,  a  vitiated  air  heater  and 
finally  to  the  head-end  assembly.  The  air  flow  can  be 
either  vented  into  the  atmosphere  or  vented  through  the 
ramjet  motor  by  the  control  of  two  pneumatically  operated 
valves.  Two  flexible  air  flow  lines  were  used  to  connect 
the  main  air  line  to  the  air  heater.  The  latter  was  mounted 
on  the  thrust  stand. 
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C.  AIR  HEATER 

To  simulate  actual  flight  condition,  air  should  be 
heated  to  the  appropriate  stagnation  temperature.  In  this 
investigation  ethylene  gas  used  in  the  air  heater.  The 
oxygen  in  the  air  which  is  consumed  during  the  heating  must 
be  replaced  before  injection  into  the  ramjet  cumbustor. 
Figure  3.2  shows  a  schematic  diagram  of  the  air  heater.  The 
air  heater  and  the  ramjet  motor  assembly  were  mounted 
together  on  the  thrust  stand. 

D.  DATA  ACQUSITION  AND  CONTROL  SYSTEM 

The  primary  instrumentation  used  in  this  investigation 
consisted  of  various  individual  pressure  transducers  and 
thermocouples  and  a  strain  gage  load  cell  for  the  thrust 
measurement.  All  transducer  outputs  were  recorded  both  with 
a  Honeywell  1508  Visicorder  and  HP-9836S  computer  using  a 
HP-3054A  AUTOMATIC  DATA  AQUISITION/ CONTROL  SYSTEM.  A  timing 
reference  signal  was  provided  to  the  analog  record  by 
feeding  a  10  Hz  signal  from  a  laboratory  signal  generator  to 
the  timing  channel  of  the  Visicorder.  During  the  test,  the 
HP-9836S  computer  system  was  used  to  record  the 
temperatures,  and  thrust,  and  calculate  the  flow  rates  of 
the  ignition  gas,  purge  gas,  air  heater  ethylene  and  oxygen, 
air  and  secondary  injection  gas.  It  was  necessary  to  record 
the  PMM  fuel  grain  weight  and  inner  diameter  before  and 
after  each  run.  The  fuel  grain  weight  was  measured  using  a 
balance.  All  thermocouples  used  were  chromel  vs  alumel 
(type  K)  with  electronic  ice  points. 

The  computer  was  used  to  both  control  the  test  sequence 
and  to  provide  data  aquisition  during  the  experiment.  Figure 
3.3  is  a  block  diagram  of  the  computer  program.  (See 
Appendix  A.  for  a  complete  listing).  Data  aquisition  for 
all  fourteen  channels  were  recorded  every  0.5  seconds. 
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I  Variable  definitions  and  nomenclature 


Transducer  calibrations 

main  air  sonic  choke  pressure 
chamber  pressure 
motor  head  pressure 

air  heater  fuel  sonic  choke  pressure 
oxygen  make-up  sonic  choke  pressure 
ignition  fuel  sonic  choke  pressure 
purge  gas  sonic  choke  pressure 
thrust  transducer 


Pre-run  input 

test  number,  date 

fuel  I.D 

fuel  dimensions  and  weight 

air  heater  fuel  type 

motor  throat  diameter 

ignition  delay  time 

ignition  time 

buring  time 

purge  time 

sonic  choke  diameters,  discharge  coefficients 

gas  contants  and  desired  flow  rates 


Flow  rates  set-up 
main  air 
heater  fuel 
make-up  oxygen 
ignition  fuel 
purge  gas 


Test  and  data  collection 

pressures  and  temperature   in  voltages 


Data  extraction 

pressures,  temperatures  and  flow  rates 


Print  results 


Figure  3.3    Block  Diagram  of  the  Computer 
Data  Aquisition/Control  Program. 
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IV.  EXPERIMENTAL  PROCEDURES  AND  TEST  CONDITIONS 

A.   TRANSDUCER  CALIBRATION 

A  dead-weight  tester  was  used  for  the  calibration  of  all 
pressure  transducers.  Initially  all  transdusers  were  checked 
for  linearity  and  the  Visicorder (analog)  display  positions 
were  labeled.  For  the  computer  data  aquisition,  voltages 
corresponding  to  atmospheric  pressure  and  the  maximum 
pressure  were  measured.  From  these  readings  a  calibration 
constant  (K)  was  determined  for  each  transducer. 


i  max, 

where     Vpmax  =   voltage   reading   at  the   applied 
maximum  pressure 

Vp0     =   voltage   reading   at   atmospheric 
pressure 

Pmax   =  maximum  applied  pressure 


The  thrust  transducer  was  calibrated  in-place  by 
appliying  known  weights  to  the  transducer  through  a  pulley 
system  attached  to  the  thrust  stand. 

B.   FLOW  RATES  SET-UP 

By  use  of  the  computer,  the  gas  flow  rates  could  be  set  up 
before  the  hot  firing.  These  were  determined  from  short 
duration  flows  using  the  equation 


30 


^^ 

where  km  =  constant  (a  function  of  R  and  |  ) 
Cd  =  discharge  coefficient 
Pt  =    total  pressure  of  the  gas  flow 
Dchoke  =  sonic  choke  diameter 
Tf  =  total  temperature  of  the  gas  flow 

C.   TEST  SEQUENCE 

1.  Thrust  Stand  Calibration 

The  thrust  stand  was  calibrated  to  determine  the 
difference  between  the  measured  thrust  and  the  actual  or 
theoretical  thrust.  The  measured  thrust  can  be  in  error 
primarily  due  to  the  variation  in  flex-line  stiffness  with 
changing  pressure  levels.  The  process  was  begun  by 
dead-weight  calibration  of  the  load  cell.  Cold  air  was  then 
flowed  through  the  motor  with  a  choked  exhaust  nozzle. 
Measured  thrust  was  recorded  and  compared  to  the  theoretical 
thrust.  The  latter  was  calculated  using  Equations  (2.11)  and 
(2.13).  This  was  repeated  for  increasing  air  flow  rates  and 
for  three  different  exhaust  nozzle  throat  diameter.  An 
additional  test  was  made  using  hot  air  (1200°  R)  and  a  12 
inches  long  HTPB  fuel  grain.  This  test  was  conducted  in 
order  to  compare  the  temperature  rise  combustion  efficiency 
determined  by  using  chamber  pressure  and  by  using  thrust. 

2.  Fuel-wall  Inj  ect ion  Tests 

For  the  fuel-wall  injection  tests,  three  injection 
location  were  used  as  discussed  above.  Each  test  lasted  for 
approximately  six  seconds,   three   seconds  without  injection 
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and  three  seconds  with  injection.  Air  injection  mass  flows 
of  2,  5  and  7  percent  of  the  main  air  flow  rate  were  used. 
For  oxygen  injection,  1.5,  2  and  5  percent  were  used.  Those 
conditions  which  resulted  in  the  greatest  increases  in 
chamber  pressure/ thrust  were  repeated  with  full  duration(30 
sec )  tests . 

3  .   Face  Inj  ect ion  Tests 

In  these  experiments  a  twelve  inches  long 
cylindrically  perforated  PMM  fuel  grain  was  used.  Air 
injection  flow  rates  of  2,  5  and  7  percent  of  the  main  air 
flow  rate  were  examined. For  oxygen  injection,  1.5,  2  and  5 
percent  were  used.  Ethylene  was  also  used  with  0.5,  1.25  and 
2.5  percent.  Testing  methods  were  identical  to  those 
discussed  above  for  fuel-wall  injection. 

4.   Inlet  Air  Swirl  Tests 

Several  inlets  were  initially  tried  which  induced 
swirl  to  the  entire  inlet  air  flow.  Swirling  the  entire  flow 
apparently  distroys  the  flameholding  recirculation  zone,  and 
resulted  in  no  ignition.  The  tube-in-hole  injector  discussed 
above  was  then  tested  using  a  twelve  inch  PMM  fuel  grain. 
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V.  RESULTS  AND  DISCUSSION 

Summaries  of  the  experimental  results  obtained  during 
the  initial  screening  tests  are  presented  in  Table  1  and 
Table  2  for  fuel-wall  and  inlet  step-face  injection 
respectively.  Results  from  the  subsequent  full  duration 
tests  are  presented  in  Table  3. 

A.   FUEL  WALL  INJECTION 

In  data  presented  in  Table  1  shows  that  air  injected 
through  the  fuel  grain  surface  was  generally  detrimental  to 
combustion.  Air  injection  rates  up  to  approximately  7  %  of 
the  inlet  air  flow  rate  had  little  if  any  effects  on 
combustion  pressure  when  introduced  downstream  of  flow 
reattachment.  Small  amounts  of  air  injected  into  the 
recirculation  zone  appeared  to  provide  small  improvements. 
The  injected  air  was  quite  cold  (approximately  500  R)  and 
therefore  could  have  had  a  quenching  effect  on  the  boundary 
layer  combustion.  Use  of  "in-flight"  high  air  should  be 
examined  in  the  future. 

Unheated  oxygen  injection  resulted  in  significant 
increases  in  combustion  pressure,  especially  when  introduced 
into  the  recirculation  zone  in  small  amounts.  Based  on  these 
screening  test  results  two  additional  tests  were  conducted 
with  test  times  of  approximately  30  seconds.  These  tests 
were  made  to  examine  the  effects  of  1  %  oxygen  injection 
into  the  recirculation  region. 

In  these  tests,  a  total  grain  length  of  fourteen  inches 
was  used;  two  inches  of  PMM  ,  then  the  injection  ring  and 
then  twelve  inches  of  PMM.  One  test (Table  3,  no.  30)  was 
made  without  injection  and  one  with  (Table  3,  no. 29). 
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Oxygen  injection  resulted  in  approximately  a  5  % 
increase  in  in  fuel  regression  rate  (when  adjusted  for  the 
different  air  flowrate)  but  a  decrease  of  21  %  in  combustion 
efficiency.  The  recirculation  zone  is  generally  fuel  rich. 
Thus,  oxygen  injection  should  increase  the  "flame-holder" 
temperature.  This  was  evident  during  the  test  from  increased 
flame  luminosity. 

It  is  not  clear  at  this  point  why  the  combustion 
efficiency  decreased  or,  for  that  matter,  why  the  precence 
of  the  ring  alone  increased  combustion  efficiency.  The 
oxygen  injection  velocity  was  approximately  30  ft/sec  which 
is  less  than  30  %  of  typical  recirculation  zone,  near-wall 
velocities.  However,  the  injected  oxygen  may  have  penetrated 
the  shear  layer. 

Further  testing  would  be  benificial,  especially  using 
HTPB,  but  at  this  point  it  does  not  appear  that  fuel  wall 
injection  of  small  amounts  of  air  or  oxygen  is  a  viable 
method  for  providing  fuel  mass  flowrate  control  since 
combustion  efficiency  was  decreased  significantly  and  the 
increase  in  regression  rate  was  small. 

B.   INLET  STEP  FACE  INJECTION 

The  data  from  the  screening  tests  with  step-face 
injection  are  presented  in  Table  2.  Unheated  air  had  no 
significant  effect  on  combustion  pressure.  Oxygen  injection 
resulted  in  increased  flame  luminosity  in  the  recirculation 
zone  but  apparently  did  not  significantly  affect  downstream 
combustion  (even  though  the  overall  equivalence  ratio  was 
approximately  0.7).  Injection  of  ethylene  resulted  in 
significant  increases  in  combustion  pressure. 

A  full  duration  test  (Table  3,  no.  28)  was  conducted 
using  1  %  of  ethylene  with  an  injection  velocity  of 
approximately  65   ft/sec.    Comparison   of  the   results  with 
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those  for  no  injection  (Table  3,  no.   27)  showed  that  only  a 
small  increase  occurred  in  fuel  regression  rate. 

Some  of  the  ethylene  apparently  escaped  the 
recirculation  region  and  burned  downstream  in  the  central 
air  or  in  the  aft  mixing  region,  resulting  in  little  effect 
on  regression  rate. 

C.  INLET  AIR  SWIRL 

Only  one  test  was  conducted  using  inlet  air  swirl  (Table 
3,  no.  31).  Approximately  43  %  of  the  air  was  injected 
axially  to  maintain  the  recirculation  zone  flame  holding 
ability.  This  could  be  reduced  to  perhaps  only  10  -  20  % , 
resulting  in  more  air  with  swirl.  The  amount  of  swirl  was 
also  intentionally  kept  small  (  5  degrees  from  the  axial 
direction)  to  determine  if  regression  rate  was  sensitive  to 
the  swirl. 

The  small  amount  of  swirl  increased  the  fuel  regression 
rate  by  15  %.  This  resulted  in  an  increase  in  equivalence 
ratio  ratio  from  0.69  to  0.79  .  Some  non-uniformity  in 
regression  rate  was  also  evident. 

These  initial  results  indicate  that  fuel  regression  rate 
is  quite  sensitive  to  inlet  air  swirl.  Further  testing  is 
necessary  using  varying  amounts  of  swirl,  but  the  technique 
appears  viable  for  in-flight  fuel  mass  flow  modulation. 

D.  THRUST  STAND  CALIBRATION 

The  results  from  the  cold-flow  tests  for  comparison  of 
measured  and  calculated  thrust  are  in  Figure  5.1.  Excellent 
agreement  was  attained.  A  least-square  fit  of  the  data 
resulted  in  the  following  equation  for  relating  measured 
thrust  to  the  "actual"  thrust 

Fa   =  1.0149  X  Fm  +  0.4  (5.1) 
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where    Fa  =   actual  thrust,    based   on  Pc   and 
throat  diameter 

Fm  =  measured  thrust 
A  test  was  then  coducted  (Table  3,  no.  32)  in  which  HTPB 
fuel  was  used  with  an  inlet  air  temperature  of  approximately 
1200  R.  This  test  was  made  in  order  to  compare  the 
combustion  efficiencies  based  on  chamber  pressure  and  based 
on  thrust.  The  two  efficiencies  were  within  5  %,  with  the 
value  based  on  Pc  being  greater.  This  difference  is  not 
large  when  it  is  realized  that  small  errors  in  determining  F 
and  Pc  from  the  analog  traces  are  squared  when  calculating  C 
for  the  efficiency  based  on  thrust  determination. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

Additional  tests  using  fuel-wall  and  step  face 
injection  of  oxygen,  gaseous  fuel  and/or  heated  air  should 
be  conducted  to  verify  the  initial  results  found  in  this 
investigation.  However,  the  initial  data  from  this 
investigation  indicates  that  gaseous  injection  is  not  a 
viable  tecnique  for  fuel  regression  rate  control. 

Fuel  regression  rate  was  found  to  be  quite  sensitive  to 
small  amounts  of  inlet  air  swirl  without  large  changes  in 
combustion  efficiency.   Additional  testing  is  recommended. 
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APPENDIX   A 
COMPUTER   PROGRAM   FOR   EXPERIMENT    CONTROL   AND    DATA   REDUCTION 

10        ISFRJTEST 

20        ! SOLID  FUEL  RAMJET  DATA  ACQUISITION  AND  DATA  REDUCTION  PROGRAM, 

30        !THIS  PROGRAM  IS  DIVIDED  INTO  FIVE  PARTS: 

40        !  (1)  VARIABLE  DEFINITIONS  AND  NOMENCLATURE 

(2)  TRANSDUCES  CALIBRATIONS 

(3)  PRE-RUN  INPUTS, FLOU  SET-UPS 

(4)  THE  TEST  SEQUENCE  AND  DATA  COLLECTION 

(5)  POST-RUN  OPERATIONS 

(1)  VARIABLE  DEFINITIONS  AND  NOMENCLATURE 


30 

60 
70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

298 

300 

310 

320 

330 

340 

350 

360 

370 

330 

390 
400 

410 
420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 


SYMBOL 

A 

Cdair 

Cdhf 

Cdho 

Cdif 

Cdo 

Cf 

Cstarth 

Cstarair 

Dairchoke 

Dates 

Dhfchoke 

Dhochoke 

Di 

Difchoke 

Do 

Doctioke 

D'th 
Dtfteff 

F 
FuelidS 

Gc 

Heaterfuel* 

Ignitionfuell 

Knair 

Kahf 

Knho 

Knif 

Km 

Lo 

Mair 

Maird 

Mhf 

rthfd 

Mho 

Mhod 

Mif 

Mifd 

Md 

Mpd 

Pa 

Pbar 

Pc 

Ph 

Phf 

Pho 

Pif 

Pd 

PuraeqasS 

Ta  ' 

TestnoS 

Thf 

Tho 

Ti 

Tid 

Tif 


DEFINITION 
ANALOG  CHANNEL  NUMBER 
DISCHARGE  COEFFICIENT, UIR  SONIC  CHOKE 
DISCHARGE  COEFFICIENT,  HEATER  FUEL  SONIC  CHOKE 
DISCHARGE  COEFFICIENT,  OXYGEN  MAKE-UP  SONIC  CHOKE 
DISCHARGE  COEFFICIENT,  IGNITION  FUEL  SONIC  CHOKE 
DISCHARGE  COEFFICIENT,  PURGE  GAS  SONIC  CHOKE 
THRUST  COEFFICIENT 
THEORETICAL  C»,  FT/SEC 
C*  FOR  AIR  FLOU  BASED  ON  Dtheff 
AIR  SONIC  CHOKE  DIAMETER 
Test  Date  Mo-Dav-Yr 
AIR  HEATER  FUEL  SONIC  CHOKE  DIAMETER 
OXYGEN  MAKE-UP  SONIC  CHOKE  DIAMETER 
MOTOR  INLET  DIAMETER,  IN. 
IGNITION  FUEL  SONIC  CHOKE  DIAMETER 
MOTOR  FUEL  PORT  DIAMETER, IN. 
PURGE  GAS  SONIC  CHOKE  DIAMETER 
MOTOR  EXHAUST  NOZZLE  THROAT  DIAMETER, IN, 
EFFECTIVE  THROAT  DIAMETER, IN. 
THRUST 

FUEL  IDENTIFICATION 
GRAVITATIONAL  CONSTANT  (32.174) 
FUEL  IDENTIFICATION 
SFRJ  IGNITION  FUEL  I  ."3 
AIR  SONIC  CHOKE  FLOy  RATE  CONSTANT 
HEATER  FUEL  SONIC  CHOKE  FLOW  RATE  CONSTANT 
OXYGEN  MAKE-UP  SONIC  CHOKE  FLOa  RATE  CONSTANT 
IGNITION  FUEL  SONIC  CHOKE  ROW  RATE  CONSTANT 
PURGE  GAS  SONIC  CHOKE  FLOW  RATE  CONSTANT 
FUEL  GRAIN  LENGTH,  IN. 
AIR  FLOW  RATE,  L3M/SEC 
DESIRED  AIR  FLOW  RATE,  Lbn/sec 
HEATER  FUEL  aQU  RATE,  LWsec 
DESIRED  HEAIER  FUEL  FLOW  RATE  Lbrt/sec 
HEATER  OXYGEN  FLOU  RATE,  Urn/sec 
DESIRED  HEATER  OXYGEN  FLOU  RATE,  Lbn/'sec 
IGNITION  FUEL  FLOW  RATE,  Lbh/sec 
DESIRED  IGNITION  FUEL  FLOU  RATE,  Lbn/sec 
PURGE  GAS  FLOW  RATE,  Lbn/sec 
DESIRED  PURGE  GAS  FLOU  RATE ,  Lbn/sec 
PRESSURE,  AIR  SONIC  CHOKE,  Psia 
BAROMETRIC  PRESSURE,  Psia 
PRESSURE,  CHAMBER,  Psia 
PRESSURE,  MOTOR  HEAD-END.  Psia 
PRESSURE,  HEATER  FUEL  SONIC  CHOKE,  Psia 
PRESSURE',  HEATER  OXYGEN  SONIC  CHOKE,  Psia 
PRESSURE,  IGNITION  FUEL  SONIC  CHOKE,  Psia 
PRESSURE,  PURGE  GAS  SONIC  CHOKE,  Psia 
PURGE  GAS  I.D 

TEMPERATURE,  AIR  SONIC  CHOKE,  R 
TEST  I.D  NO. 

TEMPERATURE,  HEATER  FUEL  SONIC  CHOKE.  R 
TEMPERATURE,  OXXYGEX  MAKE-UP  SONIC  CHOKE,  R 
TEMPERATURE,  MOTOR  AIR  INLET,  R 
TEMPERATURE,  DESIRED  MOTOR  AIR  INLET,  R 
TEMPERATURE.  IGNITION  FUEL  SONIC  CHOKE,  R 
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690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

610 

820 

830 

848 

850 

868 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

983 

990 

1000 

1010 

1020 

1330 

1040 

1053 

1060 

1070 

1080 

1090 

1100 

It  10 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310  T 

1320 

1330 

0)t560 

1340 

7) +  660 


T«a 

Tub 

T« 

TrtO 

Tp 

Utf  i 


TIME  BETWEEN  SFRJ  AIR  FLOW  AND  IGNITION,  sec 
BURN  TIME  AFTER  Tni.  sec 
IGNITION  DURATION  AFTER  T«a,  sec 
PURGE  DERATION  AFTER  Tub.  sec 
TEMPERATURE,  PURGE  GAS  SONIC  CHOKE,  R 
INITIAL  FUEL  GRAIN  WEIGHT,  Lbn 


BEEP  1000, .1 

PRTNT  USING  "9" 

PRINT  USING  '6/' 

PRINT  "  1.2  Volts  HAXIMUH  into  Acquisition  Systen' 

IHE  RECORDED  VARIABLES  (VOLTAGES)  AND  LOCATIONS  ARE: 
(NOTE:  THE  MAXIMUM  ALLOWABLE  VOLTAGE  INTO  THE  SYSTEM  IS  1.2  VOLTS) 


VARIABLE 

Pa 

Pc 

Ph 

Phf 

Pho 

Pif 

Pp 
F 

Ta 

Thf 

Tho 

Ti 

Tif 

To 


3497  DACU  SCANNER  NUMBER  0 


0 

1 

3 

4 

6 

7 

5 

9 

10 

11 


13 


ALL  FLOW  RATES  ARE  CALCUIATED  USING  THE  ONE-DIMENSIONAL,  ISENTPOPIC 
FLOiJ  EXPRESSIONS  WITH  FIXED  PROPERTIES.  SHALL  Sf'NIC  NOZZLES  HAVE 
MEASURED  DISCHARGE  COEFFICIENTS.  THE  AIR  FLOW  NOZZLE  USES  AN  ASSUMED 
DISCHARGE  COEFFICIENT  (Cd)  OF  9.97. 

rt  (LBH/SEC)=Cd*P*A*K«/T\5 

Kn  IS  THE  GAS-DEPENDENT  SONIC  CHOKE  FLOW  RATE  CONSTANT 
Kn=SQR((Gan«a*Gc/R)»(2/(GaMMa+l))A((GaHna+l)/(Ganna-l))) 

APPPOPRIATE  CONSTANTS  ARE: 

GAS    MOLECULAR  WT.   GAS  CONST.    CP    GAMMA     Kn 


AIR 

02 

CH4 

C2H4 

ARGON 


29.0 

32.0 

16.03 

28.03 

39,9 


53,3 
43,3 

96.4 

cc  i 
jj.  i 

38.7 


240 

1.40 

217 

1.45 

,5589 

593 

1,32 

.3376 

400 

1.22 

.4935 

124 

1.67 

.  uulO 

ALL  THERMOCOUPLES  USED  ARE  CHROMEL  vs  AlUMEL  (TYPE  K)  WITH 
EIECTRONIC  ICE  POINTS,  TEMPERATURE  READINGS  (VOLTAGES)  ARE 
CONVERTED  TO  DECREES  RANKINE  (R)  PER  'INDUSTRIAL  INSTRUMENTATION"  BY 
D.P,  ECKMAN  (PAGE  369).  THIS  CALCINATION  IS  PERFORMED  IN  SUBROUTINE 
Teak.  TEN  VOLTAGE  INTERVALS  ARE  USED  BETWEEN  460  AND  1460  R, 
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Ko=0       !  POST  RUN  CONTROL  FLAG 

GOTO  Transcal 

calc:    ! 

IF  Volts(. 00153  THEN  T=( (Volts*. 00068)/, ( 
IF  Volts>=. 00153  AND  Volts*. 00382  THEN  T= 


!00220)+460 
((Volts-. 00151)/. IC( 


IF  Volts>=. 00382  AND  Volts(. 00609  THEN  T=( (Volts-, 00382)/. 000022 
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1350 

1-60 

5)+B60 

1370 

?)+9M) 

1359 

3)+1069 

1390 

4)+116Q 

noo 

6) +1260 

1410 

71+1360 

1420 

1430 

1440  T« 

1450 

1460 

1475 

1480 

i  m 

1500 

1510 

1520  U 
1530  ! 

1540 
1550 
1560 
1570 
1580 
15?0 
1600 
1610 
1620 
163C 
1640 


I-  V'lt?}=. 03609  AND  Volts< 

IF  Vol ts>=. 00831  AND  Volts( 

IF  yolts>=. 01056  AND  Volt5< 

IF  Volts>=. 01285  AND  V«lts( 

IF  Volts)=. 01513  AND  Volts< 

IF  Volts>=. 01752  AND  Volts( 

IF  Volts>=. 01989  AND  WoltsC 


00331  THEN  T=< (Volts-. 00609)/. 000022 
01056  THEN  T=((Volts-. 00831)/. 000022 

01285  THEN  T=( (Volts-, 01056)/, 000022 
01518  THEN  T=< (Volts-. 01295)/, 000023 
01752  THEN  T=((Volts-, 01513)/. 000023 
01988  THEN  T=((Volts-. 01752)/, 000023 
02225  THEN  T=((Volts-. 01988)/, 000023 


RETURN 
ohiqh : 


IF  Volts), 02225  THEN  GOTO  Toohiqh 


PRINT  "TEMPERATURE  SET  AT  1460  R" 

PRINT  'I=*,I 

T=l 

BEEF  500,. 1 

WArT  .1 

KE?  "00,  ,1 

seal:  ' 

•S3)  JRW'.SE'JCER  CALIBRATIONS 

ITHERE  APE  6  PRESSURE  TRANSDUCERS  AND  ONE  LOAD  CELL  FOR  TH'IJST  WHICH 
!fi:JST  EE  CALIBRATED.  TRANSDUCER  LINEARITY  MUST  BE  VERIFIED  BEFORE  THIS 
'CALIBRATION  PROCEDURE  IS  EMPLOYED,  THE  ORDER  OF  CALIBRATION  IS  AS 
IFDLIWJS:  Pa.  Pc  Ph.  Phf,  PSo.  Pif,  Po,  F 

!IHE  FOLLOWING  Tw'u  LINES  SET  UP  722  AND  709  FOR  DATA  ACQUISITION" 
CLEAR  709        !  CLEAR  3497A  DATA  A5UISITI0N/C0NIR0L  UNIT 
CLEAR  72?        i  CLEAR  3456A  DIGITAL  VOLTMETER 
REtlOTE  709        !  MAKE  3497A  UNIT  TO  REMOTE  MODE 
OUTPUT  722;'L1R31STNZ110ST1T4QX1"  !  10  READING  PER  TRIGGER  AND  STORE 
1650  Pacal:  ! 

1660   ' >«»** *S X **K**** * ************ ******** ***<*******¥********!(******* »*»*#» 

1670  PRINT  '«  CALIBRATION  OF  Pa,  THE  AIR  SONIC  CHOKE  PRESSURE  TRANSDUCER" 
1630  U********************************************************************* 

1690  PaOcal:  " 

1700    PRINT  "  *****  ZERO      PRESSURE  *****" 

PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

DISP  "HIT  CONTINUE  WHEN  READY  TO  TAKE  ZERO  READING" 

PAUSE 


1710 
1720 
1730 
1740 


1760 
1770 
1730 
1790 
1800 
1310 


AQUIRE  CHANNEL  2 

INI TATE  SINGLE  TRIGGER 
READ  VOLTAGE 


REMOTE  709 

OUTPUT  709;"AC2" 

.'AIT  2 

OUTPUT  722;"T3" 

ENTER  722;Voa0 

PRINT  'Voa0=";Vpa0 

BEEP 

INPUT  "READING  OK?  (Y/N>",ZzS 
1820  IF  Zzt="N'  THEN  GOTO  PaOcal 
1830  Pansxcal:  ! 

1840  PRINT  USING  "g*     !  C1EAR  SCREEN 
1350  PRINT  "   *****  CALIBRATION  t***x" 
1960  PPINT  "APPLY  MAXIMUM  PRESSURE  USING  IKE  DEAD-WEIGHT  TESTER" 
1370  INPUT  "ENTER  THE  MAXIMUM  PRESSURE  IN  p?ig\Pa«ax 
1330  DISP  'HIT  CONTINUE  WHEN  READY" 
1390  PAUSE 
1900  REMOTE  709 
1910  OUTPUT  709;"AC2" 
1923  WAIT  2 
1930  OUTFUT  722; "13" 
1940  ENTER  722;Voa*ax 
1950  PRINT  ■Vcanax=";Vpa*ax)"Panax=  'jPanax 
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I960  Koa=^a«ax)/(VPi 

1970  PRINT  "Kpa=  ";Kpa 

19S0  BEEP 

1990  INPUT  -READING  OK?  <Y/N)\Zz$ 

20BO  IF  Zi$="N'  THEN  GOTO  Par.axcal 

2010  Pccal:! 

2020  ixxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 

2330  PRINT  "**CALIP.RATION  OF  Pc,  THE  SFRJ  MOTOR  CHAMBER  PRESSURE  TRANSDUCER**' 

2040  ixxx<xxx*xx*xxxx**x*xx»xxx-**xxx**xxx*x*xx**xxxxxx***xxxxxx*x*x**xx*xx**t 

2050  PcOcal:  ! 

2060  PRINT  "xxx*  ZERO  PRESSURE  •  ***■ 

2070  PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER' 

2080  DISP  'HIT  CONTINUE  WHEN  READY' 

2090  PAUSE 

2100  REMOTE  709 

2110  OUTPUT  709;'AC0' 

2120  WAIT  2 

2130  OUTPUT  722; "T3" 

2140  ENTER  722;Uoc0 

2150  PRINT  "Vpc0="jVpc0 

2160  BEEP 

2170  INPUT  "READING  OK?  <Y/N>\Zz$ 

2180  IF  Zz$='H*  THEN  GOTO  PcOcal 

2190  Pcnaxcal:  ! 

2200  PRINT  USING  '@" 

2210  PRINT  "  xxxx  CALIBRATION  **»x" 

2220  PRINT  'APPLY  THE  MAXIMUM  PRFSSURE  USING  DEAD-WEIGHT  TESTER" 

2230  INPUT  "ENTFR  THE  MAXIMUM  PRESSURE  IN  psig'.Pcnax 

2240  DISP  'HIT  CONTINUE  WHEN  READY' 

2250  PAUSE 

2260  REMOTE  709 

2270  OUTPUT  709; "ACO" 

2280  WAIT  2 

2290  OUTPUT  722;"T3" 

2300  ENTER  722;(lpc«ax 

2310  PRINT  ,Vpc«ax=";UDC«ax,"Pc«ax=,;Pciiax 

2320  Koc=Pr.nax/(VDe«ai-VpcO) 

2330  PRINT  "Kpc=  w:Kpc 

2340  BEEP 

2350  INPUT  "READING  OK?  <Y/N>\Zz$ 

2360  IF  Zi$='N"  THEN  GOTO  Pcaaxcal 

2370  Phcal:! 

2380  l*X*XXSXX*X*XX»*X*XXXXXXXSX***XXXXXXXXXX****XXX**X*X»X*«*X***«*X*X»**X-..    II 

2390  PRINT  ^CALIBRATION  OF  Ph,  THE  SFRJ  MOTOR  HEAD-END  PRESSURE  TRANSDUCER**' 

2400  IX*X»XXXXXXXXX«XXXXXXXXXXXXXXXXXXXXX*XXXXX*XXXXiiX»*<XXXXXXXXXXXXXXXXXXi-  11 

2410  PhOcal:  ! 

2420  PRINT  "xxxxZERO  PRFSSURE****' 

2430  PRINT  "INSURE  THAT  NO  PRFSSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

2440  DISP  "HIT  CONTINUE  WHEN  READY" 

2450  PAUSE. 

2460  REMOTE  709 

2470  OUTPUT  709; 'ACT 

2430  WAIT  2 

2490  OUTPUT  722; "T3" 

2500  ENTER  722;Vph0 

2510  PRINT  "Vph0=*;Uph0 

2520  EFEP 

2530  INPUT  'READING  OK?  <Y/N>'.Zz$ 

2540  IF  Zz$="N'  THEN  GOTO  PhOcal 

2550  Phnaxcal:  ! 

2560  PRINT  USING  "g" 

2570  PRINT  "****CALIBRATION****" 

2530  DISP  "APPLY  THE  MAXIMUM  PRESSURE  USING  DEAD-WEIGHT  TESTER' 

2590  INPUT  'ENTER  THE  MAXIMUM  PRESSURE  IN  PsiqVPhwax 

2600  DISP  "HIT  CONTINUE  WHEN  READY" 

2610  PAUSE 

2620  REMOTE  709 

2630  OUTPUT  709:'AC1" 
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2640  WATT  2 

2650  OUTPUT  722; "T3" 

2660  ENTER  722j0ph«ax 

2573  PRINT  'VphNax^jVphAax.'Phna^'iPt.nax 

2680  Koh=PhMx,'(Vph?iax-VphO) 

2690  PRINT  ,Koh=,;Kph 

2700  BFEP 

2710  INPUT  "READING  OK?  <Y/N>'.Zz$ 

2720  IF  Zzt='N'  THEN  GOTO  Phaaxcal 

2730  Phfcal: ! 

2740  i  «**x**»*»***»»x)(****<*«»»-«»*»s-jt-*xi{*******«-»*»»#t»*»»**»«**  «»*»««*»*!» 

2?50  PRINT  '^CALIBRATION  OF  Phf,  THE  AIR  HEATER  FUEL  PRESSURE  TRANSDtT.FR**: " 

2760  !**  **  x***Xttft*xx*K*****xx**«»*»*******  *******  ******* *********  ****h**»u 

2770  PhfOcal:  ! 

2730  PRINT  'xxxxZERO  CALIBRATION*-***" 

2790  PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

2930  PISP  "HIT  CONTINUE  WHEN  READ/" 

2910  PAUSE 

2320  REMOTE  709 

2830  OUTPUT  709;"AC3" 

2840  OUTPUT  722; *T3" 

2850  ENTER  722;0phf 0 

2360  PRINT  ,|JDhf9=,;yphfO 


2670  INPUT  "READING  OK  1     <Y/N>\Zz$ 
2830  IF  Zz$="N"  THEN  GOTO  PhfOcal 


2390  Phfaaxcal 

2900  PRINT  USING  "§" 

2910  PRINT  ac***CAL TBRATIDN***-* ■ 

2920  PRINT  "APPLY  THE  MAXIMUM  PRESSURE  USING  DEAD-WEIGHT  TESTER" 

2930  INPUT  "ENTER  THE  hAXIKUri  PRESSURE  IN  p9ig',Phfiiax 

2940  DISP  "HIT  CONTINUE  UHEN  READY' 

2950  PAUSE 

2960  REMOTE  709 

2970  OUTPUT  709;"AC3* 

29S0  OUTPUT  722; "T3" 

2990  ENTER  722j0phf«ax 

3008  PRINT  'VphfMax=":'v;phfrtax,"Phfaax=";Phfnax 

3010  Kphf=Phfttax/('v,phfnax-Vphf3) 

3020  PRINT  ,Kphf=,;Kphf 

3030  BEEP 

3040  INPUT  "READING  OK?  <Y/N>\Zz* 

3050  IF  Zz$="N"  THEN  GOTO  Phfuaxcal 

3060  Phocal:! 

3070  i*****xx********************************x***********x********************* 

3080  PRINT  "^CALIBRATION  OF  Pho,  THE  AIR  HEATER  OXYGEN  PRESSURE  TRANSDUCER** H 

3Q9fl  1  ***il***************3t*****************X*************3t******Jl************** 

3100  PhoOcal:  ! 

3110  PRINT  •*x**2ERD  CAt IBRATIGN***** 

3120  PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

3130  DISP  "HIT  CONTINUE  WHEN  READY" 

3140  PAUSE 

3150  REMOTE  709 

3160  OUTPUT  709;"AC4" 

3170  OUTPUT  722;"T3" 

3130  ENTER  722;'v'pho0 

3190  PRINT  ,VphoO=';Vphofl 

3200  INPUT  "READING  nk?  <Y/N>" ,Zz$ 

3210  IF  Zz$="N"  THEN  GOTO  PhoOcal 

3220  Phonaxcal:  ! 

3233  PRINT  USING  '&• 

3240  PRINT  "****CALIBRA7ION****" 

3250  PRINT  "APPLY  THE  MAXIMUM  PRFSSURE  USING  DEAD-WEIGHT  TESTER" 

3260  INPUT  "ENTER  THE  MAXIMUM  PRESSURE  IN  psig\Phonax 

3270  DISP  "HIT  CONTINUE  UHEN  READY" 

3280  PAUSE 

3290  REMOTE  709 

3330  OUTPUT  709;"AC4" 

3310  OUTPUT  722;"T3" 
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33?0  ENTFR  722; '>'nho«ax 

3330  PRINT  "VphoHaxs'jVphoHax  •Pho«ax=,;PhoMx 

3140  Kaho=Phc«ox/(Vpho«ax-VphoO) 

3350  PRINT  "Kpho=';Kpho 

3360  BEEP 

3370  INPUT  'READING  OK?  <Y/N>',Zz$ 

33S0  IF  1H-'H'  THEN  GOTO  Phonaxcal 

3390  Pifcal: ! 

3400  i»*xxx***'t»»»*»**«**x-»**»x*x*»x*«x***«****xx$**x»**»#**<»»*****x*****»»*x* 

3410  PRINT  "^CALIBRATION  OF  Pif,  THE  SFRI  IGNITION  FUEL  PRESSURE  TRANSDUCER** " 

3420  i**»*****»*H»#****»*»*#*#**»**###**x*»*****«**»»*****«**#******#**»**** 

3430  PifOcal:     ! 

3440  PRINT  -****ZERO  CALIBRATION****' 

3450  PRINT  'INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

3460  DISP  "HIT  CONTINUE  WHEN  READY" 

3470  PAUSE 

3480  REMOTE  709 

3490  OUTPUT  709;,AC6' 

3500  OUTPUT  722; 'T3* 

3510  ENTER  722;0pif0 

3520  PRINT  'Vpif0=";VpifO 

3530  INPUT  "READING  GIO  <Y/N)',Zz$ 

3540  IF  Zz*=*N"  THEN  GOTO  PifOcal 

3550  Pifnaxcal:  ! 

3560  PRINT  USING  "?" 

3570  PRINT  "»x#*CALIBRATION*m*" 

3580  PRINT  "APPLY  THE  MAXIMUM  PRFSSURE  USING  DEAD-UFIGHT  TESTER' 

3590  INPUT  "ENTER  THE  MAXIMUM  PRESSURE  IN  psiq'.Pifwx 

3600  DISP  "HIT  CONTINUE  WHEN  READY' 

3610  PAUSE 

3620  REMOTE  709 

3630  OUTPUT  709; 'ACS" 

3640  OUTPUT  722; "T3" 

3650  ENTER  722;Clpif«ax 

3660  PRINT  "VoifMax=";Vpifnax,'Pifnax=";Pifnax 

-  7.670  Knif=Pifnax/(Vpif«dx-VpiFB) 

3680  PRINT  "Kpif=";Kpif 

3690  BEEP 

3700  INPUT  "READING  OK?  <Y/N)\Zz$ 

3710  IF  Zz$="N"  THEN  GOTO  Pifnaical 

3720  Ppcal:  ! 

3730  i **»#**»***** x«»*if»»«»»*»**»«*xx»*»******«»*»***»x**»»«***»«*»****x*» 

3740  PRINT  "^CALIBRATION  OF  Pp,  THE  SFRJ  PURGE  CAS  PRESSURE  TRANSDUCER**" 

3750  i*-x-»x*»#xx*»**«***-**x*xft«**«*»*»«»*»«**«»**«x**»x»**«»**»»*»*»#**x*** 

3760  PoOcal:  ! 

3770  PRINT  'x***ZERO  CALIBRATION****' 

3780  PRINT  'INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

3790  DISP  "HIT  CONTINUE  WHEN  READY" 

3600  PAUSE 

3810  REMOTE  709 

3320  OUTPUT  709;'AC7" 

3330  OUTPUT  722;"T3" 

3840  ENTER  722;VppO 

3850  PRINT  'VppO^iVppO 

3860  INPUT  "READING  OK?  <Y/N)'.Zz$ 

3870  IF  Zz5=*N"  THEN  GOTO  PpOcal 

3830  Pcnaxcal:  ! 

3890  Print  using  •§• 

3900  PRINT  "k**CALIBRATION***x' 

3910  PRINT  "APPLY  THE  MAXIMUM  PRESSURE  USING  DFAD-WEIGHT  TESTFR" 

3920  INPUT  'ENTER  THE  MAXIMUM  PRESSURE  IN  psig*  Ppnax 

3930  DISP  "HIT  CONTINUE  WHEN  READY" 

3940  PAUSE 

3950  REMOTE  709 

3960  OUTPUT  709;'AC7* 

3970  OUTPUT  722; "T3" 

3960  ENTER  722;0pp«ax 

3990  PRINT  "Uopnax=,;VpciHax1,Ppnax=';Pofiax 


48 


4000  KoD=Pp«ax/(Vpp«ax-VopO) 

4010  PRINT   "Xpp=";Kpp 

4029  BEEP 

4020  INPUT   "READING  OK?  <Y/N)",Z7$ 

4040  IF  Zz$=*N"  THEN  GOTO  Pp«axcal 

4050  Fcal;  I 

40a0    '  IF  Ff=0.  THEN  GOTO  Endcal 

4070  i*********x**x  ********************  xmx********* ********************** 

4030  PRINT  "^CALIBRATION  OF  F,  THE  THRUST  STAND  LOAD  CELL**" 

4090  !H**********f*****#**********»*****XK*»H**«*****x»H**#***»*»*#***){ 

4100  FOcal:  ! 

41  IB  PRINT  "****ZERO  CALIBRATION****' 

4120  PRINT  "INSURE  THAT  THE  'ZERO'  TARE  WEIGHT  IS  ON  THE  TRAY' 

4130  DISP  "HIT  CONTINUE  WHEN  READY" 

4143  PAUSE 

4150  REMOTE  709 

4160  OUTPUT  709;"AC5" 

4170  OUTPUT  722;"T3" 

4180  ENTER  722;VfO 

4190  PRINT  "Vf0=";VfO 

4200  INPUT  "READING  OX?  (Y/N>".Zi$ 

4210  IF  Zz$="N"  THEN  GOTO  FOcal 

4228  Fnaxcal:  ! 

4230  PRINT  USING  "9" 

4240  PRINT  "****CALIERATION****" 

4250  PRINT  "APPLY  MAXIMUM  UEIGHT  TO  THE  TRAY" 

4260  INPUT  "ENTER  THE  MAXIhUM  UEIGHT  IN  LBS',F«ax 

4278  DISP  "HIT  CONTINUE  WHEN  READY" 

4280  PAUSE 

4290  REMOTE  709 

4300  OUTPUT  709;"AC5" 

4310  OUTPUT  722;"T3" 

4320  ENTER  722;(lf«ax 

4330  PRINT  ,Vfiux=";yf«ax,,F«ax=,;Fnax 

4340  Kf=Fwx/(yf«ax-yfC) 

4350  PRINT  "Kf=":Kf 

4360  BEEP 

4370  INPUT  "READING  OK?  <Y/N>\Z;$ 

43S0  IF  Z;t="N'  THEN  GOTO  Fnaxcal 

4390  Endcal:  ! 

4450  PRINT  "THIS  ENDS  THE  CALIBRATIONS  " 

4410  U-it****** *****************************************  x*********i(X'***<-»*x* 

4420  !(3)  PRE-RUN  INPUTS,  FLOW  RATE  SET-UP 

4430  i  ***** *************** ***********************iHt****x*******»**x-****»*** 

4440  !A.  PRE-RUN  INPUTS 

4450  !SM«*ff*»»*******»*»****#«******««********************************** 

4460  Inputvariables:  ! 

4470  INPUT  "DO  YOU  WANT  TO  INPUT  NEW  VARIABLES  <Y/N)",Yy* 

4430  IF  Y»»="N"  THEN  GOTO  Nochanqe 

4490  PRINT  "INPUT  THE  FOLLOUING  TEST  VARIABIES' 

4500  PRINT  "STRING  VARIABLES  MUST  BE  ENCLOSED  IN  QUOTATION  MARKS" 

4510  INPUT  "THE  TEST  IDENTIFICATION  NUHEER.A  STRING  VARIABLE", Testnot 

4520  INPUT  "TODAY'S  DATE  AS  MO, DAY. YEAR ,A  STRING  VARIABLE  ",Date$ 

4530  INPUT  "THE  BAROMETRIC  PRESSURE", Pbar 

4540  INPUT  "THE  FUEL  IDENTIFICATION, A  STRING  VARIABLE  VF«elid$ 

4550  INPUT  "THE  AIR  HEATER  FUEL  TYPE, A  STRING  VARIABIE  '.Heaterfueli 

4560  INPUT  "THE  SFRJ  IGNITION  FUEL  TYPE, A  STRING  VARIABLE  ", Iqni tionf uelt 

4570  INPUT  "THE  PURGE  GAS  TYPE, A  STRING  VARIABLE  ".Purqeqast 

4590  INPUT  "THE  INITIAL  FUEL  GRAIN  WEIGHT  (LBM)",14tfi 

4590  INPUT  "THE  INITIAL  FUEL  GRAIN  INTERNAL  DIAMETER  (IN)", Dp 

4600  INFUT  "THE  INITIAL  FUEL  GRAIN  LENGTH  (IN)',Lp 

4610  INPUT  "THE  AIR  INLET  DIAMETER  (IN)'.Dl 

4620  INPUT  "THE  MOTOR  THROAT  DIAMETER  (iN)VDth 

4630  INPUT  "THE  DESIRED  AIR  INLET  TEMPERATURE  (R)VTid 

4640  INPUT  "THE  DESIRED  IGNITION  DELAY  TIME  (SEttVTna 

4650  INPUT  "THE  DESIRED  IGNITION  TIME  (SEC)\Tni 

4660  INPUT  "THE  DESIRED  BURN  TIME  (SEO'Jnb 

4673  INPUT  "THE  DESIRED  PURGE  TIME  (SEOMno 
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46.80  INPUT  "THE  AIR  CHOKE  DISCHARGE  COEFFICIENT"  .Cdau 

4690  INPUT  "THE  AIR  HEATER  FUEL  CHOKE  DISCHARGE  COEFFICIENT" .Cdhf 

47C0  INPUT  "THE  AIR  HEATER  OXYGEN  CHOkE  DISCHARGE  COEFFICIENT" ,Cdho 

4710  INPUT  "THE  IGNITION  GAS  CHOKE  DISCHARGE  COEFFICIENT"  ,Cdif 

4720  INPUT  "THE  PURGE  GAS  CHOKE  DISCHARGE  COEFFICIENT", Cdp 

4730  INPUT  "THE  AIR  HEATER  FUEL  FLOW  RATE  CONSTANT, Knhf '  .Knhf 

4740  INPUT  "THE  IGNITION  GAS  FLOU  RATE  CONSTANT,  Knif ,Knif 

4750  INPUT  "THE  PURGE  GAS  FLDU  RATE  CONSTANT,  K«c\K*p 

4760  INPUT  "GAMMA  FOR  HEATER  FUEL  .Gan^hf  ,Ga«hahf 

4770  INPUT  "GAMMA  FOR  IGNITION  FUEL,  Ca^aif  \GaMfiaif 

4790  INPUT  "GAMMA  FOR  PURGE  GAS.  GaMaD",Gai«ap 

4790  INPUT  "GAS  CONSTANT  FOR  HEATER  FUEl  ,Rhf  (ft  ,lbf/lbn.R)\Rhf 

4800  INPUT  "GAS  CONSTANT  FOR  IGNITION  FUEL ,Rif (ft . lbf/lbn.R) ",R if 

4810  INPUT  "GAS  CONSTANT  FOR  PURGE  GAS,  Pp(ft .lbf/lh«.R)",Ro 

4820  INPUT  "THE  AIR  CHOKE  DIAMETER, Dairchoke  ( in)",Dairchoke 

4330  INPUT  "THE  HEATER  FUEL  CHOKE  DIAMPTER  .Dhtfuchokedn)"  .Dhfchoke 

4840  INPUT  'THE  IGNITION  OXYGEN  CHOKE  DIAMETER, Dhtoxchekeiin)" .Dhochoke 

4350  INPUT  "THF  IGNITION  FUEL  CHOKE  DIAMETER  ,Dianfuchok?(in)\l$ifchake 

4860  INPUT  "THE  PURGE  CHOKE  DIAMETER, Dournecfiokel in)' .Docheke 

4370  INPUT  "THE  DESIRED  MASS  FLOU  RATE  OF  AIR  ,M.urd(lb«/sec)',Maird 

4880  INPUT  "THE  DESIRED  MASS  FLOU  RATE  OF  heater  FUEL , rthf d ( lbn/sec > ' . Hhf d 

4890  INPUT  "THF  DESIRED  MASS  FLOU  RATE  OP  heater  OXYGEN, Mhod(lb*/sec)"  .Mhod 

4900  INPUT  "THE  DESIRED  MASS  FLOU  RATE  OF  IGNITION  FUEL,Mifd( lbs/sec  )'.Nif d 

4910  INPUT  "THE  DESIRED  MASS  FlOU  RATE  OF  PURGE ,Mr.d( lbs/sec)", Mod 

4920  INPUT  'THE  INJECTION  GAS , In jectqa-," , In ]ec Tq^sS 

4930       UX*XXXXXXXXXXX*X*XX*  XXXXXXX  XXXXXXXXXXXXXXXXXXXXX  XXXXXXXXXXXXXXXXX  XXXXXX 

4*40  !*   CHECK  THE  INPUT  DATA    * 

4950   IXXXXXX**X*XX*XXX«XXX-**X*XXXXXXXXXXXXX*XXXXXXX«XXXXXXXXX»XXXXXXXXX*XXXXX 

4960  Chancel:  ! 
4970  PRINT  USING  "8" 

4980  PRINT  USING  •14A,9A15X,14A,9A.5X,14A,9A,2X,14A,DDDD.DDDD";"Te5tna=,iTestno 
$,"Date=";Date$,"Fuelid$=Vuelid$,fcPbar=MiPbar 

4990  PRINT  USING  "14Ai9A,5X,14A,9A,5X,14A,9A,5X,14A,9A";  'Heaterfuel=";Heaterfoe 
1$,  "Iqniticnfuel=";Iqnitionfuel$, "Puroe  oas=" :Purneoasi 

5000  'PRINT  USING  al4A,DDDDDDD,DI;DD,5X,l4A, DDDD. DDDD, 5X,14A, DDDD, DDDD" ;"Wtfi=";W 
tfi/Dp=":Dp,'LD="]Lp 

50 1 0  PRINT  USING  M4A,DDDD,DDDD,5X,14A)Dr>DD,DDDD";"Di  =  ";Di  ,"Dth=";Dth 
5020  PRINT  USING  "  1 4A ', DDDD . DDDO , 5X ,  1 4 A  DD0D.DDDD.5X  ,14A . DDDD . CDDO ,2X ,  14A.DDDD.D 
DDD"  j  "Tm="  jTm  ,  "T«i  =  •  j  Tni ,  'Tnbs* :  Tub ,  "T«n  =  "  ;Tnp 

5030  PRINT  USING  "14A, DDDD, DDDD, 5X,14A, DDDD. DDDD";"DairchQke=";Dairchoke,"Dhtfu 
choke=';Dhfchoke 

5040  PRINT  USING  "14A, DDDD, DDDD, 5X,14A,DDDD.DDDD";"Dhtoxchoke=';Dhochoke,"D)cnf 
uchoke=*:Difchoke 

5050  PRINT  USING  ,14A)DDDD,DI;DD';*Dpurnechoke=";Dpchoke 

5060  PRINT  USING  "14A, DDDD. DDDO, 5X,14A, DDDD, DDDD, 5X ,14A,DDDD,BDDD";"Cdair=a;Crta 
ir."Cdh="jCdhf  "Cdho=*;Cdho 

507G  PRINT  USING  "14A, DDDD, DDDD, 5X,14A. DDDD, DDDD";"Cdn="  .Cdn  ,"Cdif=";Cdif 
5030  PRINT  USING  "14A, DDDD, DDDD, 5X.14A, DDDD. DDDD, $X,14A,DDPD, DDDD"; "G in»whf=" ;G 
annahf > "Gan«aif=" ;Gan«aif 

5090  PRINT  USING  "14A,DDDD.DDDDVGamiap=\Cawiap 

5100  PRINT  USING  "14A, DDDD, DDDD, 5X,14A, DDDD, DDDD,5X,14A, DDDD, DDDD"; "Rhf=";Rhf, " 
Rif=";Rif  ■Rp=';Rp 

5110  PRINT  USING  "14A, DDDD. DDDD, 5X,14A, DDDD. DDDD, 5X. HA.DDDD.BDDD" ; "K«hf =" ;K«hf 
•K«if="jK«if,"KHp=";K«p 
5120  PRINT  USING  "MA.DDDD, DDDD, 5X,14A, DDDD. DDDD";"Knho=";Knho,"Knair  =  ";Kfiai- 
5130  PRINT  USING  "14A, DDDD. DDDD, 5X,14A, DDDD. DDDD" ; "Haird=l;Haird  "Tid=";Tid 
5140  PRINT  USING  "14A, DDDD. DDDD, 5X.14A, DDDD, DDDD. 5X.14A, DDDD. DDDD.SX.HAfDDID. I) 
DDD" ;  "Mhfd="  :Mhfd .  "Mhod="  ;Mhod ,  "Mifd=" ;Mifd,  "Mpd=;;,1pd 
5150  PRINT  'INJECTION  GAS  IS",Injectaas$ 

5160   IXKXXXXXXXXX*XX-XXXXXX*XXXXXXXX*XX"XXXXXXXXXXX*XXXXXXXX*XXX»X*X*X-XXifX*«X»i  a 

5170  !x  VARIABLE  CORRECTION 

5180   !XXXXXX*XX*X*XXXXXXXXX»;SXXXXXXXXXmXXXXXXXXXXXX*XXXXX»XXXX***XX**X**X»sS-;. 

5190  PRINT  USING  "4/" 

5200  PRINT  "CHECK  AIL  VALUES" 

5210  DISP  "HIT  'CONTINUE'  TO  PROCEED  TO  CORRECTION  OR  NEXT  SUBROUTINE" 

5220  PAUSE 

5230  INPUT  "VARIABLE  OK7  <Y/N)",Zz$ 

5240  IF  Zz$="N'  THEN  GOTO  Change  variable 
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re- 
cc 


5250  IF  Zz$=T  THEN  GOTO  Nochange 

5260  Ch^naevariable: ' 

5270  PRIST  "INPUT  VARIABLE  NAHE=  CORRECTED  VA!UE" 

5230  PRINT  "ft  STRING  VARIABLE  MUST  BE  ENCLOSED  IN  QUOTATION  MARKS* 

5290  DISP  "HIT  EXECUTE  AND  THEN  CONTINUE  AFTER  CORRECTION  ■ 

5100  PAUSE 

5310  GOTO  Change! 

5320  Nochange:! 

5333  U********«t«*«jf#**#x***#»*«**»»***#K**«i{*»**»«*»*»»#»»#*««*»******»****x* 

5340  !  B.  FLOU  RATE  SET-UPS 

5350  i***«*«****»*M*»***********«*»»****»**»***M*****H****»***#*«»***#«** 

5360  PRINT  USING  'V 

5370  INPUT  "DO  YOU  WANT  TO  PRESET  THE  AIR  FLOW  RATE?U/N>\Zz$ 

5330  IF  ZiS-'N"  THEN  GOTO  P^skjp 

5390  PRINT  "SET  THE  DESIRED  VALUE  OF  Pa(psia)  USING  THE  HAND  LOADER  /PRESSURE  G 

AGE" 

5400  PRINT  USING  '3/' 

5410  PRINT  "THE  HAND  LOADER  SHOULD  BE  20  PSIG  MORE  THAN  DESIRED  PRESSURE" 

5420  Paset:  ! 

5430  PRINT  USING  "3/" 

5440  PRINT  "KANUALl Y  INITIATE  AIR  FLOW  BY  TURNING  'MAIN  AIR'  TO  'ON'  AND  PUSHI 

NG  'PR1'  ON  CONTROL  PANEL" 

5450  DISP  "HIT  CONTINUE  WHEN  READY" 

5460  PAUSE 

5470  WAIT  3 

54S0  OUTPUT  709;"AC2" 

5490  OUTPUT  722;  9T3" 

5500  ENTER  722;Vpa 

OUTPUT  70VAC9" 

OUTPUT  722;"T3° 

ENTER  722;Vla 

540  OUTPUT  7fl?;"AC0u 

5550  OUTPUT  722;"T3" 
5560  ENTER  722: Ope 
5573  OUTPUT  709 /AC3" 

5551  OUTPUT  722;'T3B 
5590  ENTER  722;Vti 
5603  OUTPUT  7C9;'AC5" 
5610  OUTPUT  722; *T3" 
5620  ENTER  722; Of 
5630  CLEAR  709 

5640  OUTPUT  709;"DC10I1"     !  DIGITALLY  CLOSE  CIRCUIT  1 

5650  WAIT  1 

OUTPUT  7Q9;"D01Q,1"     !  DIGITALLY  OPEN  CIRCUIT  1 

BEEP 

5630  PRINT  "TURN  OFF  'MAIN  AIR'" 

5690  DISP  "HIT  CONTINUE  TO  PROCEED" 

5730  PAUSF 

5710  Pa=(Vpa-VpaO)*Kpa+Pbar 

5720  Pc=(Vpc-VpcQ)*Kpc+Pbar 

5730  Volts=Vta 

5740  GQSUB  Teak          !  CONVERSION  FROM  VOLTAGE  TO  TFMPERATURE 

5750  Ta=T 

5760  Volts=Vti 

5770  GOSUB  Teak          !  CONVERSION  FROM  VOLTAGE  TO  TEMPERATURE 

5780  Ti=T 

5790  Mair=Knair*Cdair*Pa*.7B54*(DairchokeA2)/(Ta\5) 

5800  PRINT  USING  "?" 

5810  PRINT  USING  B5A>2X1DDD.DDDDD";"Mair=";Mair 

5820  PRINT  USING  "14A,d6d.DDD" ;"Mair  DESIRED=" jMaird 

5830  Raho=Mair/Maird 

5840  PRINT  USING  "20 A , D . DDD »2X ,3A , IX . DDDD . DD , 1 A , 3X ,3A , DDHD . DD . 1  A" ; "hair/DESIRE 
D  Mair=";Rano,'Ta=";Ta,'Ryh=i;Ti;,r 

5350  Pg=Pa-Pbar 

5360  PRINT  USING  "4A,DDDD.DD,4A"  ;"Pa  =';Pq;"Psia" 

5870  PRINT  USING  "4A.DDDD.DD, 1A'  "Ta  =';Ta  T 

5880  IF  Ko=l  THEN  GOTO  Dtheffcalc          !  DETERMINE  PRE-RUN  OR  POST-RUN 

5890  INPUT  "IS  AIR  FLOU  RATE  ACCURATE  ENOUGH?  (Y/N)",XxJ 
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5900  IF  Xx$="Y"  THFN  GOTO  Dtheffcalc 

5910  Par.eu=(P3*Kaird/Hair)-Pbar 

5920  PRINT  "RESET  Pa  TG';Paneu;  "Psiq" 

5930  DISP  "MIT  CONTINUE  AFTER  RCSET  GF  Pa" 

5940  PAUSE 

5950  GOTO  Paset 

5960  Dtheffcalc:          !  CAICUIATE  EFFECTIVE  THROAT  DIAMETER,  CHARACTERS 

TIC  EXHAUST  VELOCITY,  THRUST  BOTH  THEORETICAL  AND  MEASURED 

5970  Cstar  =  (GcxRair»Ti«((Ga«Maair+l)/2)M(GamiaairM)/(GaMnaair-1))/Gannaair)* 

.5 

5980  Cfair=.7396+.5283-(14.7/Pc) 

5990  Dtheff^(<(Cstar*Mair)/(.7e54*Pc*Gc))\5) 

6000  PRINT  "Cstarair(BASED  ON  Ti)=';Cstar, '(ft/sec)1 

6010  PRINT  "Cfair=,;Cfair 

6020  PRINT  "Dtheff=*;Dtheff ,"(in)' 

6030  PRINT  ■Pco=":Pc',"Psia" 

6040  PRINT  "Ti="  Ti , *  R' 

6050  F=(Vf-Vffl)*Kf 

6060  Fair=Cfair*Pc*.7854x(Dthefr2) 

6370  PRINT  "Fair(BASED  ON  Pc)  =  ";Fair,'(Lbf ) ' 

60S0  PRINT  ■Fair<MFASURED)=";F1'(Lbf)" 

6090  IF  Kc=0  THEN  GOTO  6170 

6100  INPUT  "DQ  YOU  WANT  PRINTOUT  OF  POST  RUN  DATA?(Y/N)',Yy$ 

6110  IF  Yy$="N'  THEN  GOTO  Finish 

6120  Headinoprint:  ! 

6130  PRINTER  IS  701 

6140  PRINT  USING  "6/" 

6150  PRINT  *           *xxxx*  POST  RUN  DATA  COLD  AIR,  CHOKED  *«**»*  ' 

6160  GOTO  Preprint 

6170  INPUT  "DO  YOU  WANT  PRINTOUT  OF  PRE-RUN  DATA?(Y/N>",Xx$ 

6180  IF  Xx$="Y"  THEN  GOTO  Preorint 

6190  GOTO  Skippnnt 

6200  Preprint:   ! 

6210  PRINTER  IS  701 

6220  IF  Ko=l  THEN  GOTO  6240 

6230  PRINT  "             «*»  PRE-RUN  DATA,  USING  COLD  AIR  ONLY  ANN  CHO 

KEO  FLOli  ***»■ 

6240  PRINT  "• 

6250  PRINT  USING  M2A22X,8A   ■•'TEST  DATE  ISVDate* 

6260  PRINT  'CstaraiHBASED  ON  Ti)=*;Cstar, "(ft/sec) ' 

6270  PRINT  ■Cfair=";Cfair 

6280  PRINT  'Dtheff^'jDtheff/dn)' 

6290  PRINT  "Pco=";Pc  "Psia" 

6300  PRINT  MTi=",Ti  '  R" 

6310  PRINT  'FaiMBASED  ON  Pc)=';Fair  "(Lbf )' 

6320  PRINT  'FaiMMEASlMD^'jF;  (Lbf)" 

6330  PRINT  ,Hair=';Mair,,Lbn/s5C* 

6340  PRINTER  IS  1 

6350  IF  Ko=l  THEN  GOTO  Finish 

6360  Skippnnt;  ! 

6370  DISP  "HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLOW  RATE  SET  UP' 

6380  PAUSE 

6390  Paskip:   ! 

6400  PRINT  USING  *@" 

6410  INPUT  'DO  YOU  WANT  TO  PRESET  THE  HEATER  FUEL  FLOU  RATE?  (Y/N)*,2:$ 

6420  IF  Zz$='N"  THEN  GOTO  Phfskio 

6430  Umx*mxx*»*x»*mx*«xx*xx»x»*»xxxx*x*«xxx*x*x»*«**«m#o*»xi*#*x***->.- 

6440  PRINT  "SET  THE  DESIRED  VALUE  OF  Phf  USING  THE  HAND  LOADER/PRESSURE  GAGE" 

6450  i*x*x*****«xxmxmxxxxxxxx»m*xxx«x**x*******xxxxxx*txx*x*»**m*xxxiuii 

6460  DISP  "HIT  CONTINUE  WHEN  READY' 

6478  PAUSE 

6430  Phfsets! 

6490  PRINT  USING  '§" 

6500  PRINT  'MANUALLY  TURN  ON  AIR  'HEATER  GAS'  SWITCH" 

6510  DISP  "  HIT  CONTINUE  TO  PROCEED* 

6520  PAUSE 

6530  OUTPUT  709;*AC3" 

6540  OUTPUT  722;"T3" 
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6550  ENTER  722;','phf 

fc5t0  OUTPUT  7Q9;"AC13' 

6573  OUTPUT  722; 'T3" 

b58l  ENTFR  722-,ythf 

6590  CLEAR  709 

6600  OUTPUT  709; "DC10 ,0" 

6610  PRINT  'MANUALLY  TURN  OFF  AIR  'HEATER  GAS'  SUITCH" 

6620  BEEP 

6630  DIS?   'HIT  CONTINUE  TO  PROCEED' 

6640  PAUSE 

6658  OUTPUT  709:"DO10)0" 

6663  Phf=(Vphf-6phfQ)»Kphf+Pbar 

6670  Volts=0thf 

6630  GOSUB  Teak 

669 J  Thf=T 

6760  Hhf=Kflhf*Cihf*Phf*.7B54*(Dhfchoke*2)/(Thf*.5) 

6710  PRINT  USING  T 

6720  PRINT  USING  ,4A.DD,DDDDD";'nnf=,;hhf 

6730  PRINT  USING  *12A,DD.DDDDD4,  'Mhf  DE3IRLD=';Hhfd 

6740  Ratio=rthf/Mhfd 

6753  PRINT  USING  " 1 3A , D . DDD , 2X , 4A , DDDD . Db , 1 A ' ; 'Mhf /  Mhf  DESIRED=';Ratio,'Thf=' 

,TLf  Mil                         '          III                               II                                                                             I                       I 

6763'  Pq=Phf-Fb*r 

6770  PRINT  USING  "SA.DDDD.DDD, 4A,3X,4A, DDDD. DD,1A.4A. DDDD. DD,lA';'Phf=  'iPg;"P 

Bin"  *J^t  z*  •  tj- r  ■  "0  •                        i     i      i     i               iii               ii              i    s, 

6790 '     INP'uf  "IS  HEATER  FUEL  FLOW  RATE  ACCURATE  ENOUGH?  <Y/N>\Xx$ 

6790  IF  t%%='V  THEN  GOTO  Phffin 

6638  Phfnew=(Phf*Mhfd/Hhf)-Pbar 

6810  PRINT  USING  M3A. DDDD. DD,4A'i 'RESET  Phf  TO'iPhtnewj'Psig' 

6823  DISP   'HIT  CONTINUE  AFTER  RESET  OF  Phf 

6830  PAUSE 

6340  GOTO  Phfset 

6850  Phffin;! 

6863  DISP  'HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLOW  RATE  SET  UP' 

6S70  PAUSE 

6833  Phfskip:   ! 

6898  PRINT  USING  "S1 

6900  INPUT  'DO  YOU  UANT  TO  PRESET  THE  HEATER  OXYGEN  FLOW  RA7E?<Y/N> " ,ZiS 

6910  IF  Z:$='N*  THEN  GOTO  Phoskic 

6920  ll***H*H**S*«»H«**H**H***»***«**«t*##*«*****»*******»«********l 

6930  PRINT  'SET  THE  DESIRED  VALUE  OF  Pho  USING  THE  HAND  LOADER /PRESSLRE  GAGE' 

6940  l  **»l«*«**»**»«*H**#M**«**#*****H»»»«***««tf*«**««****f  **»»** 

6950  DISP  "HIT  CONTINUE  WHEN  READY' 

6968  PAUSE 

6970  Phoset: ! 

6983  PRINT  'MANUALLY  TURN  ON  AIR  'HEATER  GAS'  SWITCH' 

6993  DISP  'HIT  CONTINUE  TO  PROCEED" 

7333  PAUSE 

7010  OUTPUT  709;'AC4' 

7020  OUTPUT  722; 'T3' 

7C3S  ENTER  722;6pho 

7040  OUTPUT  739;'ACir 

7050  OUTPUT  722;'T3" 

7063  ENTER  722;Vtho 

7373  CLEAR  709 

7083  OUTPUT  709; "DC10.0" 

7090  PRINT  'MANUALLY  TURN  OFF  AIR  'HEATER  GAS'  SWITCH" 

7138  BEEP 

7113  DISP  'HIT  CONTINUE  TO  PROOFED ■ 

7120  PAUSE 

7130  OUTPUT  7C9i "DQ10 ,0" 

7140  Pho=<Vphc-Vpho8)*Kpho+Pbar 

7153  Volts=Vtho 

7163  GOSUB  Teak 

7178  Tho=T 

7183  Hho=K*ho*Criho*Pho«.78j4*(DhochokeA2)/(Tho\5) 

7193  PRINT  USING  '?' 

7208  PRINT  USING  ,4A,DD.DDDDD';,Hho=,;Mho 
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7210  PRINT  USING  ,llA/DD.DDDDD,;"hairinjeci=,;f1ho 

7220  PRINT  USING  "18A ,DD. DDDDD" ; "hha  DESIRED=';Hnod 

7230  Ratio=H.ho/Mho<J 

7240  PRINT  USING  *33A,  DDDD . DDDD , 3X , 4A , DDDD . DD , 1 A  * ; "Mh  o/Mh  o  DE3IRED='  .Rat 

io,  ■Tho=":Tho,*R" 

72$0  Pg=Pho-Pbar 

7260  PRINT  USING  "5A ,DDDD . DD , IX ,5A ,5X, 4A ,DDDD . DD , IX ,2ft" ; "Pho=" ;Pq ; "Psiq " 

7^270  "    'INPUT  'IS  THE  HEATER  OXYGEN  FIQU  RATE  ENOUGH?  (Y/N)?',Xx$ 

72S0  IF  Xx$='Y"  THEN  GOTO  Phofin 

7290  Phonew=(Phc*Hhod/Hho)-Pbar 

7300  PRINT  USING  "1 4A , DDDD . DD , 1 X , 4A" : "RESET  Phj  TO  " ;Phonew; 'Psiq " 

7310  DI3P  'HIT  CONTINUE  AFTER  RESET  OF  Pho" 

7320  PAUSE 

7330  GOTO  Phoset 

7340  Phofin;! 

7350  DISP  'HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLQU  RATE  SET  W 

7360  PAUSE 

7370  Phoskip:   ! 

7380  PRINT  USING  '(?' 

7390  INPUT  "DO  YOU  WANT  TO  PRESET  THE  IGNITION  FUEL  FLOW  RATE?(Y/N>',Zz$ 

7400  IF  Zz$='N-  THEN  GOTO  Pifskip 

7410  imm»}u-*xxx**x»*»xx*xx*xxm»»*»*x****m»*»»x»*xxm»***»xx»»**xx«»** 

7420  PRINT  "  SET  THE  DESIRED  VALUE  OF  Pif  USING  THE  HAND  LOADER /PRESSURE  GAGE" 

7430  Uxx*x»xmm»xxx»x*mxxxx»x«»Mmo4mxm*»**x»*»**x»M»xx*»»xx****x* 

7440  DISP  'HIT  CONTINUE  WHEN  READY' 

7450  PAUSE 

7460  Pifset: ! 

7470  PRINT  'MANUALLY  TURN  ON  '  IGN.  GAS  '  SWITCH' 

7480  DISP  '  HIT  CONTINUE  WHEN  READY* 

7490  PAUSE 

7500  OUTPUT  709; 'ACS' 

7510  WAIT  1 

7520  OUTPUT  722; "T3' 

7530  ENTER  722; Up  if 

7540  OUTPUT  709;'AC12" 

7550  OUTPUT  722; 'T3' 

7560  ENTER  722;0tif 

7570  BEEP 

7590  PRINT  'MANUALLY  TURN  OFF  'IGN  GAS'  SWITCH" 

7590  DISP  'HIT  CONTINUE  TO  PROCEED" 

7600  PAUSE 

7610  Pif=('Jpif-VpifO)*Kpif+Pbar 

7620  ViltssVtif 

7630  GGSUB  Tcalc 

7640  Tif=T 

7650  Hif=K«if»Cdif*Pif»,7854*(Difchoke*2)/(Tif\5) 

7660  PRINT  USING  "?H 

7670  PRINT  USING  ■4A.DD.DDDD0" : "Mif^" ;hif 

7690  PRINT  USING  "15A,DD.DDDDD';  'hif  DESIRED=';Hifd 

7690  Ratio=Hif/Mifd 

7700  PRINT  USING  '17A,D.D&D,3X,4A,DDDD,DDD,lA';'Hif/Hif  DESIRED=';Ratio,'Tif=' 

■  Ti  +  D  " 

^710 '    Pg=Pif-Pbar 

7720  PRINT  USING  '4A)DDDD.DD,lX)4A(5X)4A,DDDD.DD)lX)lA';'Pif=';Pg;"Pisq";'Tif= 

7730  '  INPUT  'IS  THE  SFRJ  IGNITION  FUEL  FLOW  RATE  ACCURATE  ENOUGH?  <Y/N)',Xx$ 

7740  IF  Xx$='Y'  THEN  GOTO  Piffin 

7750  Pifnew=(Pif*Nifd/Mif)-Pbar 

7760  PRINT  USING  '12A>DDDD. DD. IX. 4A"; "RESET  Pif  TO" ;Pifnew;'P5iq' 

7770  DISP  'HIT  CONTINUE  AFTER  RESET  OF  Pif" 

7780  PAUSE 

7790  GOTO  Pifset 

7800  Piffin: ! 

7810  DISP  "HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLOW  RATE  SET  UP" 

7820  PAUSE 

7830  Pifskio:       ! 

7840  PRINT  USING  "0" 
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7850  INPUT  "DO  YOU  WANT  TO  PRESET  THE  PURGE  FLOU  RATE?  <Y/H>",Zi$ 

7869  IF  Z:$='N'  THEN  GOTO  Ppskip 

7S30  PRINT  'SET  THE  DESIERD  VALUE  OF  Pp  USING  THE  HANDLOA OER/ P RESSUR E  GAGE" 

7390  i  *************  t**m******s*m!u*m************'*****m********<******* 

7900  DI5P  'HIT  CONTINUE  WHEN  READY' 

7910  PAUSE 

7920  Ppser:! 

7930  CLEAR  709 

7940  OUTPUT  709;  "DCt0,l " 

7950  WAIT  2 

7960  OUTPUT  709;"AC7" 

7978  OUTPUT  722;"T3' 

7980  ENTER  722jUpo 

7990  OUTPUT  709;*feC13" 

8000  OUTPUT  722;"T3" 

8010  ENTER  722:Vtp 

8020  OUTPUT  7C9; '1)010,1' 

8030  Pp=(Vpp-Vpp0)*Kpp+Pbar 

SG40  Volts=Vtp 

8050  GOSUB  Teak 

8060  Td=T 

8070  Mp^«p*Cdp*Pp*.7354*(DpchokeA2)/(Tp\5) 

SG8I  PRINT  USING  *?' 

3090  PRINT  USING  "3A ,DD . DDDDD" ; "«p= ' ;Mp 

8100  PRINT  USING  s13A,DD.DDDI)Di;,hp  DESIRED=D;Mpd 

8110  Ratio=Hp/Hpd 

8128  PRINT  USING  ,16A,D.DD,3X)3A,DDDD.DD;lAa;,ttp/Hp  DESIRED  =" ;Ratio , "Tp=" ■ Tp 

dl3B  BEEP 

8140  Pq=Pp-Pb^r 

8150^  t   PRINT  USING  "3Af DDDD . DD, IX , 4A,3X ,4A ,BI>DD . DD,lXf 1A" ; "Pp=" jPgj "Psig" , "Tp=" 

dl60  INPUT  "IS  THE  PURGE  GAS  FLOW  RATE  ACCURATE  ENOUGH?  <Y/N>\Xx$ 

8170  IF  Xx$="Y"  THEN  GOTO  Pp.fin 

8130  Ppnew=(Pp*Hpd/Mp)-Pbar 

8190  PRINT  USING  *1 tA . DtDD . DD , 1 X ,4A ■ : "RFSET  Pp  TO" ;Ppnew; "Psiq" 

8200  DISP  "HIT  CONTINUE  AFTER  RESET  OF  Pp' 

8210  PAUSE 

8220  GOTO  Ppset 

8230  Ppfin:  ! 

8240  Ppskip:    ! 

8253  PRINT  'THIS  COMPLETES  PRE-RUN  SET-UP' 

8260  i  **********)»  ********************** ******** **************  ***xx****x*m*t*x* 

82^0  !  (4)  THIS  PORTION  OF  THE  PROGRAM  RUNS  THE  TEST  AND  COLLECTS  THE  DATA 

3230  ' ********************************** *x* ****** J*x**»x**t*x***x*x***x* ******** 

8290  ! 

8300  PRINT  USING  '?' 

8310  DISP  'SET  TIMEDATE  BY  PRESSING  K19  AND  UPDATE,  THEN  EXECUTE,  THEN  HIT  CON 

TINUE" 

8320  BEEP 

8330  PAUSE 

8340  CLEAR  709 

8350  CLEAR  ',0'D 

8360  !  THE  FOLLOWING  PROGRAMS  THE  3456A  DM 

3370  OUTPUT  722; •L1Z0L0F1 OSTDPOFL 0R3 . 1 0ST 1 1STNS01T4QX1 " 

8330  PRINT  'THE  DESIRED  RUN  TIMES  AND  DELAYS  NCU  BE  INPUT" 

3390  INPUT  'ENTER  TnaJniJ/ibJup  SEPARATED  BY  COMMA. THEN  HIT  CCNTIUE",Tna,Trti 

.Tub, Tup 

8400  Tshut=Tna+Tni+T«b+Tnp+l     !  ONE  MORE  SECOND  FOR  TAIL  PART  READING 

8410  T4=T«a+Tfu 

8420  T5=T4+T*b 

8430  T6=T5+Tnp 

8440  OPTION  BASE  1 

8450  DIM  In(600,14)           !  STORAGE  DEFINITION 

8460  J=l 

8470  A=l                   !  CONTROL  FLAGS 

8430  A1=0 
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8490   A2=0 

8500   A3=0 

8510   A4=0 

8520   VI  4=0 

8530   PRINT  USING  '8' 

8540   DIS?  "TO  INITIATE  RUH,  PUSH  'PRI-AIR  SWITCH'  ■ 

8550  Monitor:    !  WAIT  UNTIL  INITIATE 

8560         OUTPUT  709;"AC14* 

8570         OUTPUT  722;"T3U 

8530         ENTER  722;  VI 4 

8590         IF  U14>=.5  THEN  GOTO  Startrun   !  IF  INITIATE  RUN,  IT  UILL  GIVE 

0.8  VOLTS 
360Q         GOTO  Monitor 
8610  Startrun:  ! 
8620        T0=TIMEDATE 
8630       ON  CYCLE  .5  GDSUB  Data 
8640        ON  DEI  AY  Tshut  GOTO  Shutdown 
8650  Tinewait:  !  UAIT  0,5  SECOND 
3660        A=A+1 
8670        GOTO  Tinewait 
8680  Data:  !   READ  ALL  DATA  EVERY  0.5  SECONDS 
8690     OUTPUT  709; "AC0AFGAL13' 
8730    FOR  1=1  TO  14 
8710     OUTPUT  722;'T3" 
8720     ENTER  722  USING  "# ,K * ; In C J , I > 
8730     OUTPUT  709; "AS" 
8740     NEXT  1 
8750     J=JM 
8760     T1=TIMFDATE 
8770     IF  Al=l  THEN  GOTO  Tnabypass 
8780     IF  DROliND(Tl-T0 ,3) >=Tna  THEN 
8790      OUTPUT  709j*DC10,2" 
8330       OUTPUT  709; "DC10,3" 
8810      A1=l 
8820     END  IF 
8830  Taabvpass:  ! 
8340     tl=TIMFDATE 
8850     IF  A2=l  THEN  GOTO  T4bypass 
8860      IF  DRCUND(T1-T0.3)>=T4  THEN 
8370      OUTPUT  709 v -DO10 ,2* 
3380       OUTPUT  709; '0010,3' 
8390       A2=l 
8900      END  IF 
8910  T4bypass:  ! 
8920     T1=TIMEDATE 
893C     IF  A3=l  THEN  GOTO  T5bypass 
8940      IF  DROUND(Tl-T0,3))=T5  THEN 
8950       OUTPUT  709;'DC1Q,1* 
8960       A3=l 
8970      END  IF 
8930  TSbypass:  ! 
8990     T1=T1MEDATE 
9000     IF  A4=l  THEN  GOTO  Tfcbypass 
9010      IF  DROUNOCT1-T0 ,3) >=T6  THEN 
9020       OUTPUT  709;,DO10,l" 
9030       OUTPUT  709;  "0010,0" 
9040       A4=l 
9050     END  IF 
9060  TSbypass:  ! 

9070        GOTO  9080   !»  THIS  NUMBER  UILL  BE  CHANGED  WHEN  RENUMBER, 
9080   RETURN 
9090  Shutdown:  ! 
9100        OFF  CYCLE 
9110         Jnax=J-t 
9120        PRINT  USING  '9* 
9130        PRINT  '  TEST  COMPLETE,  TURN  OFF  'MAIN-AIR',  AND  TURN  OFF  'HFATER 

CASES' " 
9140         BEEP 
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9150  !*#a»*x«»*****»K********K****«K#*«»**#*«**»*»***********»»*#**x****»****** 

9160  !*  (5)  POST-RUM  OPERATION 

9180        DISP  °HIT  CONTINUE  TO  PROCEED  TO  DATA  REDUCTION' 

9190        PAUSE 

9200        OUTPUT  709;'DO10,0" 

9210        PRINTER  IS  701 

9220  PRINT  USING  '3/' 

9230  PRINT  '  ****  PRE-RUN  INPUT   ***»  ' 

9240  PRINT  USING  '2/' 

9250  PRINT  USING  ,14A,9A,5X,14A>9A,5X,14A,9A,2X,14A>Dr.DD.DDDD";,,Testno=";Tesinc 

$/»ate=,;Date$,,Fuelid$=4;F»elid$,iPbar=,;Pbar 

92i0  PRINT  USING.  "14A!9A)5X)14A,9A,5X,14Ar9A,5X)14A,9A";,Heaterfuel=";He3te-f^ 

1$, ■Iqnitionfuel=*iIanitionfuel$f "Purge  gas=";PyroeQas$ 

9273  "PRINT  USING  ' 1 4A , DDDD . DDDD , 5X , 1 4A , DDDD . DDDD , 5X , 1 4A , DDDD . DDDC ' ; "Utf i=" j U  r i 

,,DD='jDo.'Lp=,;Lp 

9230  PRINT  USING  "1 4A ,DDDD . DDDD , 5X , 1 4A , DDDD . DDDD" : "Di=" ; Di  ,"Dth=";Dth 

92<?0  PRINT  USING  'HA. DDDD, DDDD, 5X, HA, DDDD, DDDD, 5X.HA, DDDD. DDDL .2X,1 4A,  11 .jj  D 

DDD" ; aT«a=' ;T«a,'T«i=" }Tm , "Tubs' ;Tnb , "Tf«=" ;Tm 

9300    PRINT  USING  'HA, DDDD, DDDD, 5X,HA. DDDD. DDDD':  "Dairchake=VDaarchoke)'Dnt  1 1, 

chBke=\Dhf  choice 

9310    PRINT  USING  "14A, DDDD. DDDD, 5X,14A, DDDD, DDDD" j "Dhtoxchoke=" jDhochoke, "Diqnf 

ucheke="jDif  choke 

9320    PRINT  USING  '14A,DDDD.DDDDVDcurqechoke=";Dnchoke 

9330    PRINT  USING  "14A, DDDD. DDDD, 5X,14A;DDDD,DCDD15X/14A, DDDD, DDDD"; "Crtaip=";Cda 

ir,"Cdh=";Cdhf  "Cdho=';Cdho 

9340  PRINT  USING  'HA, DDDD, DDDD, 5X. HA, DDDD, DDDD\'Cdp=" :Cdn,'Cdif=';Cdif 

9350  PRINT  USING  "HA, DDDD. DDDD. 5X, HA, DDDD. DDDD, 5X,14A, DDDD, DDDD";"Gartftahf=" ;G 

a««ahf  ,  °Gap.naif=*;Gamaif 

9360  PRINT  USING  •14A,DDDD.DDDD,;,GaMa0=,;Caflfnp 

9370  PRINT  USING  "HA, DDDD. DDDD, 5X, HA, DDDD. DDDD, 5X,14A, DDDD, DDDD" ; "Rhf=" jRhf , ■ 

Rif-*'Rif  " d ->  =  •  •  p o 

9330  'PRINT  kit   "HA, DDDD. DDDD, 5X, HA, DDDD. DDDD, 5X, HA, DDDD, DDDD" ; "Krthf =" ;K«hf 

, aKnif=";K»if , "K«p=" jKnp 

9393  PRINT  USING  "HA, DDDD, DDDD, 5X, HA, DDDD, DDDD";'K«ho=';K«ho,"K«air=";Knair 

9400  PRINT  USING  "HA, DDDD. DDDD, 5X, HA, DDDD, DDDD " j "Maipd=4 ;Ma ird  "Tid=" ;Tid 

9410  PRINT  USING  "HA, DDDD. DDDD, 5X. HA, DDDD. DDDD, 5X, HA, DDDD, DDDD2X, HA, DDDD. D 

DDD":"«hfd="jMhfd  "hhod=";Hhod  "Hifd="  :flifd,  "Mpd=fc;Kpd 

9420  PRINT  'INJECTION  GAS  IS  \Injectq8S$ 

9430  PRINT  USING  "3/" 

9440  PRINT  '  ****  DATA  EXTRACTED   **** 

9450  PRINT  "" 

9460  PRINT  "FLOW  RATES  IN  Lbn/sec,  PRESSURE  IN  Psia,  TEMPERATURE  IN  R" 

9473  PRINT  " 

9430  i  HEADING  PRINT 

9490  PRINT  USING  •9A,1X.6A,3X,6A,4X.6A,4X,6A,4X,6A.4X,6A,5X16A,4X16A,4X,1CA':,T 

ir,8(sec)\"  Mair"  "    Hhf1,'  Hho"'    flif1,"    Mp1,      ti",*F(Lbf )',"  Pc","  Ph","ta" 

9500    FOR  J=l  TO  Jnax 

9510    FOR  1=1  TO  14 

9520       IF  1=1  THEN 

9530         Vpc=In(J,I) 

95^0         Pc=(Vpc-0pc0)*Kpc+Pbar 

9550      END  IF 

9560      IF  1=2  THEN 

9570         Vph=In(J,I) 

9580         Ph=(Voh-to>hO)«Kph+Pbar 

9590      END  IF  ' 

9600      IF  1=3  THEN 

9610  Voa=InlJ.I) 

9620  Pa=CJpa-0pa0)»Kpa+Pbar 

9630   END  IF 

9640   IF  1=4  AND  Ht=l.  THEN 

9650  Uphf=In(J,I) 

9660  Phf=(Vphf-VphfO)*Kphf+Pbar 

9670   END  IF 

9630   IF  1=5  AND  Ht=l.  THEN 

9690  Vpho=In(J,I) 

9738  Pho=(ypho-UDhoOXKpho+Pbar 
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9710  END  IF 

9720  IF  1=6  AND  Ff=l  THEN 

9730  yf=In(J,I) 

9740  F=(Vf-VfO)*Kf 

9750  END  IF 

9760  IF  1=6  AND  Ff=fl  THEN 

9770  F=0 

9730  END  IF 

9790  IF  1=7  THEN 

9300  Vnif=In(J,I) 

9810  Pif=(Vpif-Vp  if 0)*Kp if 

9320  IF  Pif <=(Ph+20)  THEN 

9330  Pif=0. 

9840  END  IF 

9850  END  IF 

9860  IF  1=8  THEN 

9870  Vpo=In(J.I) 

9380  Po=('Jpp-0pp0)»Kpp+Pbar 

9890  IF  Pp<=(Ph+20)  THEN 

9900  Pp=0, 

9910  END  IF 

9920  END  IF 

9930  IF  1=9  THEN 

9940  Volts=In(J,I) 

9950  GOSUB  Teak 

9960  Ti=T 

9970  END  IF 

9980  IF  1=10  THEN 

9990  V«lts=In(J,I) 

10000  GOSUB  Teak 

10016  T^=T 

10020  END  IF 

10030  IF  1=11  AND  Ht=l,  THEN 

10040  Volts=In(J,I) 

10050  GOSUB  Teak 

10060  Thf=T 

10070  END  IF 

10030  IF  1=12  AND  Ht=l.  THEN 

10090  Volts=In(J,I) 

10100  GOSUB  Teak. 

10110  Tho=T 

10120  END  IF 

10138  IF  1=13  THEN 

10140  Volts=In(J,I) 

10150  GOSUB  Teak 

10160  Tif=T 

10176  END  IF 

10180  IF  1=14  THEN 

10190  Volts=In(J,I) 

10200  GOSUB  Teak 

10210  Tp=T 

10220  END  IF 

10230  NEXT  I 

10240  T=J/2 

10250  Hair=K*air»Cdair*Pa»,7854*(DairchokeA2)/(Ta\5) 

10260  Hhf=XMhf*Cdhf*Phf*.7854*(Dhfchoke*2)/(Thf\5) 

10270  Kho=Knhc*Cdho*Pho*.7354*(DhochokeA2)/(Tho\5) 

13280  Hif=KKif*Cdif*Pif*.7S54*(Difchoke'2)/(Tif\5) 

10290  Hp=Krip*Cdp*Pp*.7854*(Dc.chokeA2)/(Tp\5) 

10300  PRINT  USING  '2X,DD.DD,3X,DD.DDD,4(3X,D.DDDDD) ;5(3X,DDDD.DD),;T;M=ir,hn? 

hoJif.Mp.TiJ.Pc.Ph.Ta 

18310  VEXT  ) 

10320  PRINTER  IS  1 

18330  PRINT  USING  '&' 

10340  INPUT  'DO  YOU  WANT  TO  HAKE  POST-RUN  AIR  CALCULATION?  <Y/N)\Yy$ 

10350  IF  Yy$='Y"  THEN  Ko=l 

10360  IF  Yy$='N'  THEN  GOTO  Finish 

10370  GOTO  Paset  !  TO  CALCULATE  THE  EFFECTIVE  THROAT  DIAMETER 


10380  Finish:   ! 

10390     PRINT  'DATA  OUTPUT  IS  COMPLETE" 

10400    DISP  'SECURE  TEST  CELL  !!!' 

10410    LOCAL   70?  !  RETURN  LOCAL  MODE 

10420  END 
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